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ABSTRACT 


The open-end resonator in suspended substrate transmission line was 
analyzed in terms of two open-end discontinuities on a line segment. A 
procedure for computing the resonant frequency and constant of propa- 
gation was demonstrated using a full wave analysis. The fringing 
capacitance was computed using an equivalent length extension model 
and a transmission line circuit model. The characteristic equation was 
derived using Galerkin’s method applied in the Fourier transform do- 
main. The calculation has been carried out in two increasing orders of 
approximation, and the results compared. The discontinuity capacitance 
at the resonator open ends was calculated for a range of line dimensions 
and substrate dielectric constants. The dispersive propagation constant 


of suspended substrate line was also calculated by the Galerkin method. 
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I. INTRODUCTION 


The analysis of planar transmission lines in the Fourier transform 
domain or spectral domain is superior to many numerical methods in the 
spatial domain. The analysis in the Fourier transform domain was used 
earlier by Yamashita and Mittra [Ref. 4] for computation of the charac- 
teristic impedance and phase velocity of a microstrip line based on a 
quasi-TEM approximation. A variational method was used in the Fourier 
transform domain to calculate the line capacitance from the assumed 
charge density. This is a low-frequency approximation neglecting longi- 
tudinal electric and magnetic fields supported by the strip transmission 
line having two dielectric media in its cross section, or inhomogeneous 
line. 

As the operating frequency is increased, dispersion characteristics of 
the inhomogeneous line become important for precise designs. This re- 
quirement has led to the full wave analysis of microstrip lines, represented 
by the work of Denlinger [Ref. 15], who solved the integral equations 
using a Fourier transform technique. The solution by his method, how- 
ever, strongly depends on the assi’.ned current distributions on the strip 
in process of solution. To avoid this difficulty and permit systematic im- 
provement of the solution for the current components to a desired degree 
of accuracy, a new method was presented by Itoh and Mittra and com- 
monly called the Spectral Domain Approach ( SDA ). In SDA, 
Galerkin’s method is used to yield a homogeneous system of equations 
to determine the propagation constant and the characteristic frequencies 


from which the equivalent circuit is derived. 














In each of these methods the Fourier transform is taken along the 
coordinate axes in the plane of the strip. By virtue of the Fourier trans- 
form domain analysis and Galerkin’s method, SDA has several features: 

e Easy formulation in the form of a pair of algebraic equations. 


® Variational nature in determination of the propagation constant. 


The Spectral Domain Approach is applicable to the following struc- 
tures: 


® Most planar transmission lines such as micvostrip, finline, and 
coplanar waveguide in multilayer configurations. 


¢ Both open ard enclosed structures. 
e Slow-wave lines with lossy dielectric materials. 


® Resonators of planar configurations. 


An efficient formulation of the SDA is achieved in the present work 
by enclosing the suspended substrate in a metallic shield enclosure. In 
this way the Fourier transform becomes the finite Fourier transform, 
calculated as a summation. Singutarities of the integrand can be avoided, 
and the truncated summations have acceptable convergence behavior. 
As applied here, the analysis assumes loss-free conductors and dielectric 
media, and also assumes infinitesimal] thickness for the strip conductors. 
The resonant frequencies are obtained bv numerically solving the char- 
acteristic equation, and equivalent open-circuit capacitances are obtained 
using the end-effect and the transmission-line models. The details of the 
analysis method will appear in Chapter II of this thesis. 

In Chapter III, discussions are presented concerning the accuracy of 
the solution for the resonant frequencies, the propagation constants, and 
the end-effect at the open end of the microstrip structure by calculation 


of the fringing capacitance. These solutions for various orders of ap- 








proximation are compared, based on the assumptions for various current 


distributions in x and z - directions. 





" 
Il. SPECTRAL DOMAIN APPROACH FOR SHIELDED 
SUSPENDED SUBSTRATE RESONATOR 
A. INTRODUCTION 
The shielded suspended substrate resonator to be analyzed is shown 
in Figure 1. A rectangular strip of width w and length | is placed on the 
suspended substrate which is separated from the ground plane by an air 
gap. The sides and the top of the structure are enclosed with metallic 
shielding walls. Thus the entire structure is considered to be the sus- 
pended microstrip resonator located in enclosed partially filled 
waveguide. It is assumed that the thickness of the strip is negligible and 
that all the media and conductors are lossless. The shielding waveguide 
has dimensions 3.2 mm X 1.575 mm X 100.0 mm, which corresponds to 
a system in actual practice. For simplicity, the strip is assumed to be 
symmetrically located. 
The operating frequency is chosen to be below cutoff of the shielding 
waveguide partially filled with substrate material, to avoid interaction 
with waveguide-mode resonances. From the formula in [Ref. 10] the 


cutoff frequency for the structure in Figure | is shown in Table 1. 


Table 1. THE CUTOFF FREQUENCY OF SUSPENDED SUBSTRATE LINE 


Dielectric constant (€,) Frequency cutoff ( GHZ ) 
2.2 ( duroid ) 44.7574 
9.8 ( alumina = AZO, ) 43.3333 

















Figure 1. Shielded suspended substrate resonator 


B. METHOD OF SOLUTION FOR RESONANT FREQUENCY 
1. Field Equations 


The analysis of field equations on strip transmission lines of the 
type considered here has been carried out by Uwano and Itoh [Ref. 2]. 
These authors treated the problem by applying Galerkin’s method of 
moments in the Fourier transform domain. Their work is summarized 


in a set of Green’s function equations as follows 
2,7), + Lz), = £, (la) 


Zyl, + ZyJ, = Ey (15) 





Where the tilde over the symbols indicates the Fourier transformed 
quantities. Here, J, and J, are the transformed z and x components of 
the currents on the strip conductor. E, and E, are the components of 
electric field tangential to the substrate surface. Z. ; Zn ; Z. ; Ze are the 


Green’s impedance functions, and are defined as the following : 


= 1 pat 2s 
wg 5) [p Zot a Z, J (1c) 
a” + B 
~~ a ~ ~ 
ee ee (1d) 
a” + B 
xz ™ Sx (le) 
= 1 2 25 
Looe ee Z.+B Z,J (If) 
a” + Bp 


y2Cts + Py3Cty 
CtoCty + Cty Ct37,9/ Vy) + Cty Ctyy,3/Vy1 a Vy3/¥y2 


WN? 


i ¥29Cty + 7,3Ct 
VopVzyOty Cty + 92172301 Cts + 92272301 Ct3 + 922 


Ct; = coth yh (1i) 
Cty = coth yt (1y) 
Ct, = coth y3d (1k) 
yp =a + B= kj (1) 





k? = wo we) (1m) 
?] 

ry = Ys (17) 
” 

Y21 = 2) (lo) 

Jy = JOE) (Ip) 

Z) = JWR (1q) 

£) = EE, (Ir) 

Hy = Hy (Is) 


Where subscripts / = 1,2,3 refer to the corresponding regions 1,2,3, w is 
the operating frequency, and ¢, and yu, are the free-space permittivity and 
permeability, respectively. In this thesis, the Fourier transform is carried 
out in the bounded region interior to the enclosed resonator. In order 
to examine the layered structures in Figure | with added bounding verti- 
cal electric and magnetic walls at x = + 1.6mm , z=+50mm , the finite 
Fourier transform should be used [Ref. 20]. There we have the two- 


dimensional finite Fourier transform pairs 


ie 1.6 50 re 
Sn ¥s Bry) -| fx, y, ee HP" dxdz 2a) 
—1.69—50 

















| a eae 
fee. 92)= Zo DVS Fan ds Bde elm tym (2) 


Xn Bin 


Where o, and B,, are the discrete transform variables for the dominant 
mode and are defined by, a, =(n— 1/2)n/a for E, even, —H, odd(in x) 
modes, and f,, =(m — 1/2)n/b for E, even, —H, odd(in z) modes. Note 
that b is the half length of the resonator and is defined to be 5 times the 
length of strip, and a is the half width of the resonator. From (2b) we 


have 





2 a : ’ 
IXY; z) = us 2 >. J (Oar J Berens elPm (2c) 
(2ab) : 
n Bm 
Equations (la,b) contain four unknowns Zz : ye , E, , and E, . Two 
unknowns E, and E, , however, can be eliminated by applying Galerkin’s 
method in the spectral domain. The first step is to expand the unknowns 


J, and J, in terms of known basis functions i and J,, 


N2 

Tn Brn) = oy On Br) (3a) 
j=l 
Ni 

Fe Bp) =)» delet n Bo 3b) 
k=] 


Where c, and d, are unknown coefficients. The basis functions must be 


chosen to approximate the true but unknown distributions on the strip. 








Siace the current is nonzero only on the strip, therefore each basis func- 
tion must be chosen so that it is nonzero only on the strip. After substi- 
tuting (3) into (la,b) and taking inner products of the resulting equations 
with the basis functions Jy and J for different values of i, this process 


yields the equations, 


NI N2 
Kinde + C2) =0 ,i=1,2,3,.,N2. (4a) 
j=l k=] 
Ni N2 
> «Pg +) K?%)q=0 ,1=1,2,3,...NI. (4b) 


Where from definition of inner products associated with the Fourier 
transform defined by (2), equations (4) lead to a set of homogeneous 


equations forming the solution matrix for the c¢ and d, . The matrix el- 


ements are 
aee aFeS 
L,I 7 5 i 
Kj = 9 » ,, Jj On BiZz (On Bm )I2A%y B m) (94) 
N=— OO m= —o90 


Ke = dy, >: > ue On Brn) Zex(On Bmx Gn Bn) (5d) 





Ki) = ay y2 Tek n BrdZx10n Brolin Pm) (Se) 


N=— Op — 
+00 to 
2,2 = ia = s 
KY a dy Y » Ji Fag fe ae On Bm xn Bm) (5d) 
N=—COm=—o0 


The right hand side of (1a,b) can be eliminated in Galerkin’s process via 
the application of Parseval’s relation, because the current JX) ° JX) and 
the field components E(x,D +T), E(x,D + T) are nonzero in the com- 
plementary regions of x. Note that the solution matrix is symmetrical. 
In order to have non-trivial values for c, , d , the determinant of the 
matrix must be zero. From this we have the solution for resonant fre- 
quency. The accuracy of the solution can be systematically improved by 
increasing the number of basis function (NI + N2) and by solving a larger 
size matrix equation. However, if the first few basis functions are chosen 
sO as to approximate the actual unknown current distribution reasonably 
well, the necessary size of the matrix can be held small for a given accu- 
racy of the solution, resulting in numerical efficiency. Moreover, the 
numbers of basis functions for Jz) and J,(z) should be equal [Ref. 18]. 
Hence the choice of basis functions is important from the numerical point 
of view. In actual computations for the dominant mode [Ref. 1], Je and 


if have been chosen to be product functions of the form: 
Ji On Bm) = J (%n)J2(Bin) (6a) 


10 











Tg§On Brn) = Fx IABrn mt (66) 
Note that (6a,b) are the Fourier transforms of 
J,4(x,2) = Jy(x)J9(z) 

Je (%Z) = J3(x)J4(z) 


Where the functional forms of J, , J, , J; , J, , are given in Figure 2. 
The coordinate forms of the current distributions shown in Figure 2 are 


the following: 


AW=s—T1+1SPI (7a) 
Jz) = + eo ae (7b) 
Jx(x) = + sin( >) (7c) 
IQ (7d) 





Taking the Fourier transform of (7) we get : 


2sin(a,w) ss 3 


ne oe 2sin(a,w) é 2{1 — cos(a,,w)} 


[ cos(a,w) — aor (8a) 


(a,w)? 


Fx Brn) = = ae (8b) 


— (Bi 














ee 
os erraat 


i riz 
Jfny=— f 
itd) 7 CONS ) 





Figure 2. Forms of assumed current distributions. 





ss 2r sin(«,,W) 

J3(0,) = a at 
(a,w) — 7 

IAB, _ cos(B) __sin(B)) (8d) 


Bl (Bn? 


Corresponding to the structure in Figure 1, w in (7) and (8) is the half 


width of the strip and / is the half length of the strip. 
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2. Procedure of solution for resonant frequency 


The field equations (1a,b) 


~ ~~ w@ 


Zi), + Lr Jy = E, 


2y7I, + Zy,J, = Ey 


were derived assuming the existence of J, , and any solution derived from 
(1a,b) must satisfy the equations (1a) and (1b) simultaneously. Since there 
are four unknowns, ¢ , d, , E, ; E. in (la,b), they can not be solved in 
general. We introduce Galerkin’s procedure by taking inner products of 


(la) with J, and likewise (1b) with Te . Then, apply the (3a,b), 


i) Tied + ae) Jaded = ) Bi, (9a) 
a, B a, B a, B 

g) Teas or ae) JZ xi » Edi (96) 
a, B a, B a, B 


By the virtue of Parseval’s theorem, the right sides of (9) are zero, so we 


have 
g) dee + de) JaZach =0 (10a) 
a, B a, B 
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og) TParke a ak) JetZ Tei a0 (105) 
a, B a, B 


In second-order approximation, equations (10) form the matrix solution 
with dimensions [2 x 2] . The matrix is symmetrical. In order to have 


nontrivial values for c , d, , the determinant of the matrix must be zero. 


a. First order approximation 
For the lowest-order approximation, here called first-order, we 
assume the current only in z-direction (J,= 0,i= 1,7 = 1), and choose the 
current distribution J,, = J,(x) x J,(z). Since we have current only in the 
z-direction ( FF + 0, uk =0Q ), then from (10) there remains only the 


equation: 


1) JaBeeka = 0 (11) 
a, B 


By using trial values of frequency, the solution f, ( resonant frequency ) 


is found as the value best satisfying (11) 


b. Second order approximation 
For the second order solution we assume current exists in both 


directions ( z and x ) and choose the current distribution, 
Jn Pm) = J (On) 2(B in) 


Fey 6p, Brn) = Fx(0,)J4 Brn) 








Since we have current in the z and x-directions ( J, ¥ 0, J, # 0), then we 


have the matrix solution 


a1) JB ov a1) Jade =0 (12a) 
a, B a, B 

i ane - 4) Jad x1 = 0 (125) 
a, B a, B 


In this case we find f, as the frequency for which the determinant of the 


coefficients vanishes. 


C. METHOD OF SOLUTION FOR PROPAGATION CONSTANT 

The increasing use of microstrip lines at microwave frequencies has 
recently created considerable interest in the study of dispersion charac- 
teristics of these lines. Until quite recently, much of the work on the 
microstrip line was based on a TEM analysis. This analysis is employed 
to calculate the static capacitance of the structure from which the char- 
acteristic impedance and the propagation wavenumber are subsequently 
derived. However, this analysis, which is necessarily approximate, is inef- 
fective for estimating the dispersion properties of the line at higher fre- 


quencies [Ref. 14]. 


1. Field equation 
The analysis technique for computing dispersion characteristics of 
shielded transmission line of the type considered here has been carried 


out by Mittra and Itoh [Ref. 3]. The formulation of the problem applied 


15 





Galerkin’s method of moments in the Fourier transform domain and is 


summarized in a set of Green’s function equations like those used above. 


Zy (Ops Bn Io(Apn) + Zo Gys Bo Seren) = Ely) (13a) 
Zens Br Ie(n) + Zul tas Brn I ely) = Eq(at,) (136) 


: j Te 
Where f is the unknown propagation constant, and a, =(n— 2): 
Coefficients Z,, Z,, Z,,,Z,, are the Green’s impedance functions and are 


defined as shown in the previous section, and 


1.6 
J at,) = | J (ne 1" dx (14a) 
—1.6 
~ 1.6 ~ . 
J(a,) = | Jae dx (146) 
—1.6 


J, and J are the transforms of strip currents J, and J,. Notice that E, and 
E. are unknown since the electric fields E,(x,D + T) and E,(x,D + T) are 
unknown for w <|x|< a, though thev are zero on the strip. The field 
equations in (13a,b) are similar to the field equations in (la,b), but with 
current distributions assumed uniform in the z-direction and propayation 


constant B is now the unknown. 











2. Procedure of solution for propagation constant 
The first step is to expand the unknowns J, and J, in terms of 


known basis functions J, and J, 


N2 

F(ct,) =) gia, a, (15a) 
j=l 

T(0t,) = y dy Jep(t,) (156) 


Where o and d, are unknown coefficients. The basis functions a and 
is must be chosen such that their inverse Fourier transforms are nonzero 
only on the strip |x|<w. After substituting (15) into (13) and taking 
inner products with the basis functions ¥e and J for different values of 


i. This process yields the equation, 


NI N2 
> Khia KUAg =O ,i=1,2,3,.,N2. (16a) 
jel k=l 
Ni N2 

Kyo, + » KPa, =0 ,i=1,2,3,.,N1. (165) 
jel k=1 


Where from definition of inner products associated with the Fourier 


transform defined by (14), the matrix elements are 


17 








+00 aq 


KNB) = 6; » Fel )Z zn BIA Cn) (17a) 
n=—oo 
+c0o 
P= dh.) Fle ZadenPmualen (176) 
n=—0o 
+oo 
KNB) = 6 > Teh p) Zed Gn, Bn IA%p (17¢) 
n=—co 
+00 
KBB) = de) TehodEaalOnPd eal) (17d 
n=~00 


The right hand side of (13) can be eliminated in Galerkin’s process via 
the application of Parseval’s relation. 

Now the simultaneous equations (16) are solved for the propa- 
gation constant f at each frequency w by setting the determinant of the 


coefficient matrix equal to zero. 


a. First order approximation. 
We choose of matrix size : N] =0, N2=1. The matrix problem 


reduces to : 


a1) JalodZa, B)J,y(a) = 0 (18) 
a, B 








b. Second order approximation. 


We choose matrix of size : NlI=1, N2=1. The matrix to be 


solved is 

C1 Y Fler B)Jy(a) + ay Y Tle@Z B)J,-1(a) = 0 (19a) 
a, p a, p 

C1 Y Jambo B)J, (a) + dy ) Fang a, B)J_1(a) =0 (196) 
a, B a, B 


Typical output values for the moment-method calculation of 
propagation constant f for various frequencies are shown in Table 3 and 


Table 4 on page 24, Figure 6 on page 29 and Figure 7 on page 30. 
D. METHOD OF SOLUTION FOR FRINGING CAPACITANCE 


In this work, the resonant frequency data of the suspended substrate 
resonators have been used for estimating the end effect at the open ends 
of suspended substrate lines using the full-wave theory. In order to get the 
value of fringing capacitance, we use two methods, the equivalent length 
extension model and the transmission line circuit model. 


The descriptions of each method are as follows : 
1. Equivalent length extension model 
From the dispersion relation the guide wavelength /, = 2 =. is de- 
rived at the resonant frequency f/, of the microstrip resonator of length /. 
Consider the open-circuited resonator whose structure is identical to 
Figure | except that the length is /' instead of /. The length /’ is deter- 


g 


, 
mined from the resonant condition of the open circuited line /’ =—-. 


o~ 
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The hypothetical extension of the suspended substrate line, which ac- 
counts for the end effect, is given by A/=/' —/. In terms of TEM model 
for transmission line characteristic impedance is defined by L 


(Henry/meter) and C, (Farad/meter), and we have equations: 


Z, = /LIC, (20a) 
B=a/LC, (206) 


From equations (20a) and (20b), the capacitance C, is defined 





C,= (21) 


where C, is capacitance/meter of line. The value of fringing capacitance 
AC is 


p pie) 


ZW 





AC=AIC, =(I'~0 


In this work we calculate the end effect for a single open end, as shown 
in Figure 3 ( for fh ) using the computer program in Appendix D. 
2. Transmission line circuit model 
The value of fringing capacitance can be calculated by using the 
approach of the transmission circuit model. The input impedance seen 
from the center of the strip resonator is calculated using formula as the 
following [Ref. 12}: 


Z,+JjZ., tan(B) 
ig = : Z,*0.0 (23) 
Z, + JZ, tan(B) 
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MICROSTRIP RESONATOR 





Figure 3. End effect on substrate resonator 


Where 


ee 
ae jow,Ac 





and AC is a fringing capacitance, / is the half of strip length, f is the 
propagation constant, obtained as explain in the previous section, and 
Z, is the characteristic impedance of strip. The formula to be used to 
calculate characteristic impedance Z, is the following [Ref. 17}: 


For0<W<a 





Z, = OV + RI {6G/W+ 1+ 4(G/W) }] 
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Where, 
= —1.7866 — 0.2035(7/G) + 0.4750(2.4/G) 


R = 1.0835 + 0.1007(7/G) — 0.09457(2.4/G) 


and c is the free-space velocity of light. This formula has been confirmed 
by variational calculation [Ref. 19], and Z, is assumed here to vary only 
slowly with frequency. The computer program for this calculation is 
shown in Appendix D. In Table 2 is shown the values of Z,, the char- 


acteristic impedance for different value of €, and the width of strip . 


Table 2, THE CHARACTERISTIC IMPEDANCE OF SUSPENDED 
SUBSTRATE LINE 




















Dielectric Impedance first order 
constant 


(B), Z, 
(é,) 


cc 
a es 2 


Impedance second 
order (f), Z, 












The true value of fringing capacitance AC can be founded by using 
trial values in (23) and looking for the minimum value of the input 


impedance Z,,. 
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I. DISCUSSION OF RFSONATOR COMPUTATION 
A. INTRODUCTION 


The calculation of the suspended stripline resonator in this work used 
numerical computations in a scale model of millimeter wave integrated 
circuits. The numerical results presented in this paper have been gener- 
ated on a main frame IBM 370 using FORTRAN 77 and compiler VS 2 
FORTRAN, since this compiler is 10 times faster than WATFOR 77. 
The solution for resonant frequencies and propagation constant was made 
by taking the absolute value of the determinant matrix. This approach 
leads to the minimum value for the system determinant rather than 
crossing the axis of determinant = 0. All output data together are shown 
in Appendix A, while the complete data values are shown in Table 3 for 
é, = 2.2 and Table 4 for e, = 9.8. The computation results for fringing 
capacitance using the transmission line circuit model agree with the 
computation result using the length extension model, as is shown by the 
very small difference between results of these two methods. 

For each of these methods we can see that as the width of strip in- 
creases the resonant frequency of a resonator of given length increases as 
shown in Figure 4 and Figure 5. The complex dispersion behavior of this 
line system is shown by the decrease of B with increasing the width of 
strip. Also the value of fringing capacitance goes up with width. This 
situation may be due to the fact that the end region will store more ac- 
cumulated charge with wider lines of the strip. There will be currents 


flowing in the end region, corresponding to the extra charge [Ref. 6]. 
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Table 3. THE DATA RESULTS FOR ¢, = 2.2 















Trans- 
mission 


Method | Width | Reso- 
nant 
fre- 
quency, 










nf 
(GHZ) (Farad) 







Table 4. THE DATA RESULTS FOR «, = 9.8 














Method | Width | Reso- 
nant 
fre- 
quency, 


Trans- 
mission 












iA (Farad) 
(GHZ) 
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B. RESONANT FREQUENCY 
For the first order approximation the modified summation it’ terms 


of a and B to (11) is represented as 


+0008©=— +00 


» » Gln» Brn) = 0 (24) 


n=—O0M=—00 


Where G(a,,8,,) is the product form of basis functions and Green’s 
impedance function. The computer programs for first and second order 
approximation for J, is shown in Appendix B. These programs calculate 
the value of the determinant with various values of strip width in order 
to get the resonant frequency, with dielectric constant for substrate : 


é, = 2.2 for duroid, and ¢, = 9.8 for alumina. The solutions for first order 


are found : 

+30 +3000 

yy », J7\(Oq» BrndZz7(%qs Bin) =0 (25) 
n=—30m=—3000 


In the x-direction we take the sum from from n=-30 to n= +30 and in 
the z-direction from m=-3000 to m= +3000. Since a higher limit does 
not change the result as shown in Table 5. 

In a similar way, we may compute the second order approximation 
with the same summation limits for a, and for 8, All data outputs for 


resonant frequency appear in Appendix A. 
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Figure 4. Resonant frequency versus width of strip for 2, = 2.2 


Table 5. THE COMPARISON FOR RESONANT FREQUENCY IN FIRST 
ORDER APPROXIMATION 


Resonant Frequency 
(GHZ) 

n=-30 to m = -3000 to oe 13.023998 

n= + 30 m= + 3000 

n=-40 to m = -4000 to 2.2 13.023498 

n= +40 m= + 4000 
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Figure 5. Resonant frequency versus width of strip for «, = 9.8 
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C. PROPAGATION CONSTANT 
The program to be used for computation of the propagation constant 
is shown in Appendix C. An initial approximation is taken as the fol- 
lowing [Ref. 7]: 
i | 


a hae (26) 


ve 





Where c is the free-space velocity of light. To find the exact value of 
propagation constant, we sweep the value of B from 1% to 200% of (26) 
in order to bring the value of the matrix determinant in (18) to zero. The 
propagation constant can be identified by the occurrence of a minimum 
value of the matrix determinant. In this case we take the summation 
limits from n=-40 to n= +40, since a higher limit does not change the 


result as shown in Table 6. 


Table 6. THE COMPARISON FOR PROPAGATION CONSTANT IN FIRST 
ORDER APPROXIMATION 


B ) 


[mead vo wera [a2 [ro [sao 
[n=-00 1 n= +100 [32 [vo [ser 






( See in Appendix A for the complete results ) 
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PROPAGATION CONSTANT 
300.0300.4 300.8 301.2 301.6 302.0 


12.802.83 12.86 12.89 12.92 12.95 12.98 13.01 13.04 13.07 13.10 
FREQUENCY (GHZ) 


Figure 6. Resonant frequency versus propagation constant for ¢, = 2.2 
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LEGEND 


PROPAGATION CONSTANT 
292.§B2.98 293.46 293.94 294.42 294.90 





9.309.37 944 9.51 9.58 9.65 9.72 9.79 9.86 9.93 10.00 
FREQUENCY (GHZ) 


Figure 7. Resonant frequency versus propagation constant fore, = 9.8 
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D. PREDICTION OF FRINGING CAPACITANCE VALUE 


The value of fringing capacitance can be found by solving equation 
(22) for the equivalent length extension model and solving equation (23) 
for for transmission line circuit model. The computer program used to 
calculate the fringing capacitance are shown in Appendix D. From the 
results in Table 3 and Table 4 we see that these two methods agree since 


there is only a small difference in their results. 
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Fringing capacitance versus width of strip for «, = 2.2 
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FRINGING CAPACITANCE *10” 


23.0 24.8 


0.70 0.76 0.82 0.88 0.94 1.00 1.06 112 118 1.24 1.30 
THE WIDTH OF STRIP (MM) 


Figure 9. Fringing capacitance versus width of strip for ¢, = 9.8 











IV. CONCLUSION 


From a computational point of view, the calculation using Green’s 
function is very effective to solve the stripline resonator problem. Actu- 
ally the computed result for the resonant frequency depends on the choice 
of the definition of the current distributions in the x-direction and the 
z-direction. In the Table 3, Table 4, Figure 5 and Figure 6 it is seen that 
the resonant frequencies increases as the width of strip increases. We do 
not see this phenomenon in a resonator without suspended substrate, as 
shown in [Ref. 1]. 

The values of the fringing capacitance found here appear to be rea- 
sonable: the values are of the order of 10-'* Farad. The error between 


the two methods of calculation is very small, providing a confirmation 


of the length extension model. 





APPENDIX A. OUTPUT DATA 


RESONANT FREQUENCY 
1. First order approximation 


For name w = 0.711 mm and ¢, 


For w 


For w 


For w 


FREQUENCY (GHZ) 


12. 909532500000 
12. 909547800000 
12. 909563100000 
12. 909578300000 
12. 909593600000 


= 1.0 mm and «cp, = 2.2 


FREQUENCY (GHZ) 


13. 023622500000 
13. 023684500000 
13. 023746500000 
13. 023808500000 
13. 023870500000 


= 1.25 mm and ¢ 


FREQUENCY (GHZ) 


13. 090996700000 


13. 091095900000 
13. 091195100000 
13. 091294300000 
13. 091393500000 
13. 091492700000 
13. 091591800000 


= 0.711 mm and c, = 9.8 


FREQUENCY (GHZ) 


9. 406851964275 
9. 406858630475 
9. 406865296675 
9. 406871962875 
9. 406878629075 


= 2.2 


= 9.9 
DETERMINANT'S VALUES 


0. 0195708908 
0. 0134981312 
0. 0063938200 
0. 0007518493 
0. 0081374571 


DETERMINANT'S VALUES 


0. 0549158342 
0. 0306886621 
0. 0076220483 
0.0151751861 
0. 0382801555 


DETERMINANT'S VALUES 


- 0922133327 
- 0595099851 
. 0294705480 
. 0022094622 
- 0342780240 
. 0658574104 
. 0986036658 


ooooo0o09o 


DETERMINANT'S VALUES 


0. 0055711373 
0. 0026056429 
0. 0003598495 
0. 0033253400 
0. 0062908285 
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For w 


For w 


2. Second order approximation 


For w 


For w 


For w 


= 1.0 mm and « 


FREQUENCY (GHZ) 


- 708079097500 
- 708178995000 
- 708278892500 
- 708378730000 
- 708478687500 
- 708578585000 
- 708678482500 


woo Oww wow 


= 1.25 mm and ¢, 
FREQUENCY (GHZ) 


- 891719662500 
- 891720994500 
- 891722326500 
- 891723658500 
- 891724990500 
- 891726322500 
- 891727654500 


wow wow 0 Oo 


= 9.8 


= 9.8 


DETERMINANT'S VALUES 


. 1101123587 
- 0734216453 
. 0367312736 
. 0000412434 
- 0366484453 
. 0733377923 
. 1100267979 


ooooo0o90 


DETERMINANT'S VALUES 


. 0013436906 
. 0009188024 
. 0004939143 
. 0000690263 
. 0003558617 
. 0007807497 
. 0012056375 


ooooo0o°o 


= 0.711 mm and fr = 2.2 


FREQUENCY (GHZ) 


12. 858237475000 
12. 858239430000 
12. 858241385000 
12. 858243340000 


= 1.0 om and ¢, 
FREQUENCY (GHZ) 


12. 966809460000 
12. 966809956000 
12. 966810452000 
12. 966810948000 
12. 966811444000 


= 1.25 mm and ¢ 
FREQUENCY( GHZ) 


13. 033089704000 
13. 033102202000 
13. 033114700000 


= 2.2 


= 2.2 


THE VALUES OF DETERMINANT 





134. 7679778482 
45. 1968072590 
44. 3743231685 
133. 9454133199 


THE VALUES OF DETERMINANT 


21. 9744171435 
12. 5409965103 
3. 1075768752 
6. 3258416755 
15. 7592591259 


THE VALUES OF DETERMINANT 





95. 3389230655 
36. 4297524373 
168. 1980546933 
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For w 


For w 


For w 


= 0.711 mm and cp = 9.8 


FREQUENCY( GHZ) 


9. 354219440000 
9. 354244230000 
9. 354269030000 
9. 354293820000 


= 1.0 mm and ¢, 
FREQUENCY( GHZ) 


9. 624875000000 
9. 625000000000 
9. 625125000000 
9. 625250000000 
9. 625375000000 


= 1.25 mm and c, = 9.8 
FREQUENCY (GHZ) 


9. 780300000000 
9. 780360000000 
9. 780420000000 
9. 780480000000 





= 9. 


8 





THE VALUES OF DETERMINANT 





838. 0000000000 
369. 8750000000 
96. 0000000000 
605. 0000000000 


THE VALUES OF DETERMINANT 


2076, 6096794134 
1067. 1692908816 

57. 7672394809 
951. 5964768065 
1960. 9218599920 


THE VALUES OF DETERMINANT 


359. 7268520610 
83. 7447038466 
192. 2324952445 
468. 2047453591 
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B. PROPAGATION CONSTANT 
1. First order approximation 


For w 


For w 


For w 


For w 


For w 


= 0.711 mm and ¢, 
FREQUENCY (GHZ) 


12. 909578300 
12. 909578300 
12. 909578300 
12. 909578300 
12. 909578300 


= 1.0 mm and ¢, 
FREQUENCY (GHZ) 








13. 023746500 
13. 023746500 
13. 023746500 
13. 023746500 
13. 023746500 


= 1.25 mm and ec, 
}REQUENCY (GHZ) 


13. 091294300 
13, 091294300 
13. 091294300 
13. 091294300 
13. 091294300 


= 0.711 mm and e;, 
FREQUENCY (GHZ) 


9. 406865000 
9. 406865000 
9. 406865000 
9. 406865000 
9. 406865000 


= 1.0 mm and ¢, 
FREQUENCY (GHZ) 








. 708379000 
. 708379000 
. 708379000 
. 708379000 
. 708379000 


www wo 





= 2.2 
BETA (RAD/METER) 








301. 793213000 
301. 796631000 
301. 800049000 
301. 803467000 
301. 806885000 


= 2.2 
BETA (RAD/METER) 


300. 933559227 
300. 968062043 
301. 002564859 
301. 037067675 
301. 071570490 


= 2.2 
BETA (RAD/METER) 








300. 448130204 
300. 482811969 
300. 517493734 
300. 552175499 
300. 586857264 


= 9.8 
BETA (RAD/METER) 


294. 290520504 
294. 427868289 
294. 565216074 
294. 702563859 
294, 839911643 


= 9.8 
BETA (RAD/METER) 


————_— 





293. 233772058 
293. 375522190 
293. 517272321 
293. 659022453 
293. 800772584 
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DETERMINANT 





0. 006551314 
0. 003314178 
0. 000016721 
0. 003176749 
0. 006300699 


DETERMINANT 








0. 043851649 
0. 016535467 
0. 010783630 
0. 038105642 
0. 065430568 


DETERMINANT 








0. 047616359 
0. 023635390 
0. 000348146 
0. 024334248 
0. 048322918 


DETERMINANT 
0. 218513360 
0. 122327343 
0. 026102604 
0. 070160849 
0. 166463009 


DETERMINANT 


0. 208247936 
0. 123411162 
0. 038539359 
0. 046367465 
0. 131309303 








9. 708379000 293. 942522716 0. 216286146 


@ 


« For w = 1.25 mm and cp = 9.8 


FREQUENCY (GHZ) BETA (RAD/METER) DETERMINANT 
9. 891724000 292. 705621320 0. 146004764 
9. 891724000 292. 850048436 0. 069357733 
9. 891724000 292. 994475552 0. 007321405 
9. 891724000 293. 138902668 0. 084032642 
9. 891724000 293. 283329783 0. 160775971 
2. Second order approximation 

*® For w = 0.711 mm and ec, = 2.2 
FREQUENCY (GHZ) BETA (RAD/METER) DETERMINANT 
12. 858241000 301. 914079040 264. 8257310433 
12, 858241000 301. 916633866 188. 2549021972 
12. 858241000 301. 919188693 111. 6834688237 
12, 858241000 301. 921743520 35. 1114309231 
12. 858241000 301. 924298346 41. 4612114998 
12. 858241000 301. 926853173 118. 0344584496 
12. 858241000 301. 929407999 194. 6083099229 
12. 858241000 301. 931962826 271. 1827659171 

® For w = 1.0 mm and ec, = 2.2 


FREQUENCY (GHZ) BETA (RAD/METER) DETERMINANT 
12. 966810000 
12. 966810000 
12. 966810000 
12. 966810000 


12. 966810000 


47. 3499570198 
25. 5508907438 
3. 7517095664 
18. 0475865116 
39. 8469974898 


301. 189557419 
301. 191275018 
301. 192992616 
301. 194710215 
301. 196427814 


*® For w = 1.25 mm and ep = 2.2 
FREQUENCY (GHZ) BETA (RAD/METER) DETERMINANT 
27. 3862465221 
15. 0664900067 
2. 7466685987 
9.5732177036 
21. 8931688989 


13. 033102000 
13. 033102000 
13. 033102000 
13. 033102000 
13. 033102000 


300. 776834261 
300. 778560641 
300. 780287021 
300. 782013401 
300. 783739781 


= 9.8 
BETA (RAD/METER) 


« For w = 0.711 mm and ec, 


FREQUENCY (GHZ) DETERMINANT 





402. 7507924038 
192. 4514179721 


9. 354269000 
9. 354269000 


294. 739297051 
294. 762060358 








For w 


For w 


9. 354269000 
9. 354269000 
9. 354269000 


= 1.0 mm and « 
FREQUENCY (GHZ) 








9. 625125000 
9. 625125000 
9. 625125000 
9, 625125000 
9. 625125000 


= 1.25 mm and e¢ 
FREQUENCY (GHZ) 








. 780360000 
. 780360000 
. 780360000 
. 780360000 
. 780360000 
. 780360000 
. 780360000 


wowuwuvnvwv ww 


294. 
294. 
294. 


= 9.8 


BETA (RAD/METER) 








293. 
293. 
293. 
294. 
294. 


= 9.8 


BETA (RAD/METER) 











784823665 
807586971 
830350278 


904601953 
951446805 
998291657 
045136509 
091981361 


. 532445612 
- 561005835 
. 589566057 
. 618126279 


646686502 


- 675246724 
. 703806947 
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17. 8618972464 
228. 1891527630 
438. 5303480889 


DETERMINANT 





311. 5213476632 
120. 8819158085 
69. 7832174509 
260. 4740500634 
451. 1905799777 


DETERMINANT 
204. 6350496775 
136. 5694649559 
68. 4983401777 
0. 4216756325 
67. 6605283907 
135. 7482716025 
203. 8415537136 








FRINGING CAPACITANCE USING LENGTH EXTENSION MODEL 


1. 


For 


For 


For 


For 


For 


For 


W 


Ww 


First order approximation 





= 0.711 mm and cp = 2.2 
RESONANT PROPAGATION 
FREQUENCY CONSTANT 
12. 9095780 301. 800050 
= 1.0 mm and cp = 2.2 
RESONANT PROPAGATION 
FREQUENCY CONSTANT 
13. 0237470 301. 002570 
= 1.25 mm and c, = 2.2 
RESONANT PROPAGATION 
FREQUENCY CONSTANT 
13. 0912940 300. 517490 
= 0.711 mm and «, = 9.8 
RESONANT PROPAGATION 
FREQUENCY CONSTANT 
9. 4068650 294. 565220 
= 1.0 mm and ¢c, = 9.8 
RESONANT PROPAGATION 
FREQUENCY CONSTANT 
9. 7083790 293. 517270 
= 1.25 mm and cp = 9.8 
RESONANT PROPAGATION 
FREQUENCY CONSTANT 
9. 8917240 292. 994480 
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FRINGING CAPACITANCE 
( FARAAD ) 








0. 00000000000000805039 


FRINGING CAPACITANCE 
( FARAAD ) 








0. 00000000000001015215 


FRINGING CAPACITANCE 
( FARAAD ) 








0. 0000000000000120243 


FRINGING CAPACITANCE 
( FARAAD ) 








0. 000000000000023461 


FRINGING CAPACITANCE 
( FARAAD ) 


0. 000000000000027951 


FRINGING CAPACITANCE 
( FARAAD ) 


0. 00000000000003185787 








2. Second order approximation 


For w = 0.711 mm and ep = 2.2 


For 


For 


For 


For 


For 


W 





























RESONANT PROPAGATION 
FREQUENCY CONSTANT 
12. 8582410 301. 921740 
= 1.0 mm and ¢c, = 2.2 
RESONANT PROPAGATION 
FREQUENCY CONSTANT 
12. 9668100 301. 192990 
= 1.25 mm and c, = 2.2 
RESONANT PROPAGATION 
FREQUENCY CONSTANT 
13. 0331020 300. 783740 
= 0.711 mm and c, = 9.8 
RESONANT PROPAGATION 
FREQUENCY CONSTANT 
9. 3542690 294. 784800 
= 1.0 mm and ¢, = 9.8 
RESONANT PROPAGATION 
FREQUENCY CONSTANT 
9. 6251250 293. 998290 
= 1.25 mm and c, = 9.8 
RESONANT PROPAGATION 
FREQUENCY CONSTANT 
9. 7803600 293. 618130 
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FRINGING CAPACITANCE 
( FARAAD ) 








0. 0000000000000080382 


FRINGING CAPACITANCE 
( FARAAD ) 








0. 000000000000010099657 


FRINGING CAPACITANCE 
( FARAAD ) 








0. 0000000000000119056853476236 


FRINGING CAPACITANCE 
( FARAAD ) 








0. 0000000000000234773 


FRINGING CAPACITANCE 
( FARAAD ) 








0. 0000000000000278197027432552 


FRINGING CAPACITANCE 
( FARAAD ) 








0. 0000000000000316946 


D. FRINGING CAPACITANCE USING TRANSMISSION LINE CIRCUIT 
MODEL 
1. First order approximation 


For w 


For w 


For w 


For w 


For w 


= 0.711 mm and ec, 
FREQUENCY (GHZ) 


12. 9095783 
12. 9095783 
12. 9095783 
12. 909578? 
12. 9095783 
12. 9095783 
12. 9095783 


= 1.0 mm and ¢, = 
FREQUENCY (GHZ) 


13. 0237470 
13. 0237470 
13. 0237470 
13. 0237470 
13. 0237470 
13. 0237470 
13. 0237470 


= 1.25 mm and « 
FREQUENCY (GHZ) 


13. 0912940 
13. 0912940 
13. 0912940 
13. 0912940 
13. 0912940 


= 0.711 mm and ec, 
FREQUENCY (GHZ) 


. 4068650 
- 4068650 
. 4068650 
. 4068650 
. 4068650 
- 4068650 
. 4068650 


wowwowwww 


= 2.2 





CAPACITANCE (FARAAD) 


oooooco°o 


2. 
CAPACITANCE (FARAAD) 


ooooo°”joo 


. 00600000000000800 
. 00000000000000802 
. 00000000000000804 
. 00000000000000806 
. 00000000000000808 
. 00000000000000810 
. 00000000000000812 


” 
& 





. 0000000000000101667 
. 0000000000000101668 
. 0000000000000101668 
. 0000000000000101668 
. 0000000000000101668 
. 0000000000000101668 
. 0000000000000101669 


= 2.2 


CAPACITANCE (FARAAD) 











IN?SEDANCE 


ooooo 30 


oooo°co°o 


. 043897341449 
. 029423400665 
. 014949733818 
. 000476340899 
. 013996778098 
. 028469625184 
. 042942194364 


IMPEDANCE 


. 000020988378 
. 000010770731 
. 000000553083 
. 000009664563 
. 000019882210 
. 000030099857 
. 000040317503 


IMPEDANCE 


0. 000000000000012043 
0. 000000000000012043 
0. 000000000000012043 
0. 00000.000000012043 
0. 000000000000012043 


0. 0000123459405 
0. 0000045583584 
0. 0000032292234 
0. 0000110168053 
0. 0000188043870 


= 9.8 


CAPACITANCE (FARAAD) 








oooo0o0ceo°o 


. 00000000000002353658 
. 00000000000002353658 
. 00000000000002353658 
. 00000000000002353659 
. 00000000000002353659 
. 00000000000002353660 
. 00000000000002353660 


= 1.0 mm and c, = 9.8 


IMPEDANCE 


oo ooo°c°cjo 


000004112742 
. 000002943926 
. 000001775109 
. 000000606293 
. 000000562522 
. 000001731339 
. 000002900155 





For w 


FREQUENCY (GHZ) 








. 7083790 
. 7083790 
. 7083790 
. 7083790 
- 7083790 
. 7083790 


wowwowwo wo wo 


= 1.25 mm and « 
FREQUENCY (GHZ) 








9. 8917240 
9. 8917240 
9. 8917240 
9. 8917240 
9. 8917240 


CAPACITANCE (FARAAD. 








ooooo0o$9o 


. 000000000000¢ 28040 
. 000000600000028045 
. 000000000000028050 
. 000000000000028055 
. 004000000000028060 
. 000000000000028065 


= 9.8 
CAPACITANCE (FARAAD) 








0. 
0. 
0. 
0. 
0. 


000000000000031975 
000000000000031976 
000000000000031977 
000000000000031978 
000000000000031979 


2. Second order approximation 


For w 


For w 


For w 


= 0.711 mm and ¢, 


FREQUENCY (GHZ) 


12. 8582410 
12. 8582410 
12. 8582410 
12. 8582410 
12. 8582410 


—_— 





= 1.0 mm and ec, = 2 


“as 


FREQUENCY (GHZ) 








. 9668100 
. 9668100 
12. 9668100 
. 9668100 
. 9668100 
. 96681090 


= 1.25 mm and ¢, 


FREQUENCY (GHZ) 








. 0331020 
. 0331020 
. 0331020 
13. 0331020 
13. 0331020 
13. 0331020 
. 0331020 


= 2.2 
CAPACITANCE (FARAAD) 


0. 
0. 
0. 








000000000000008044 
000000000900008046 
000000000000008048 


0. 00000000000C 308050 


0. 


000000000000008052 


2 


CAPACITANCE (FARAAD) 





. 000000000000010111 
. 000000000000010112 
. 000000000060010113 
. 000000C000e 7010114 
. 000000000000010115 
. 000000000000010116 


ooooo°cjo 





= 2.2 
CAPACITANCE (FARAAD) 


ooooo°oceo 








. 00000000000001192344 
. 00000000000001192345 
. 00000000000001192346 
. 00000000000001192347 
. 00000000000001192348 
. 00000000000001192349 
. 00000500000071192350 


IMPEDANCE 





. 002439262912 
. 001337532326 
. 000227305924 
. 000877916293 
. 001983634327 
. 003089348177 


qQqo0o ooo 


IMPEDANCE 
0. 000409506598 
0. 000233952816 
0. 000058399157 
0. 000117154378 
0. 000292707792 





IMPEDANCE 





0. 003037869907 
0. 001608804997 
0. 000179742742 
0. 001249316856 
0. 002678373799 


IMPEDANCE 
. 001425888209 
. 000922250353 
. 000418612916 
. 000085024103 
. 000588660705 
.001092296E£88 


ooo0o0o°o 


IMPEDANCE 


. 000010349874 
. 000006513849 
. 000002677825 
. 000001158198 
. 000004994223 
. 000008830247 
. 00001.666272 


qooooo0oo 





For w 


For w 


For w 


= 0.711 mm and cy, = 9.8 

FREQUENCY (GHZ) CAPACITANCE (FARAAD) 
9. 3542690 0. 0000000000000235503 
9. 3542690 0. 0000000000000235506 
9. 3542690 0. 0000000000000235509 
9. 3542690 0. 0000000000000235512 
9. 3542690 0. 0000000000000235515 
= 1.0 mm and c, = 9.8 

FREQUENCY (GHZ) CAPACITANCE (FARAAD) 
9. 6251250 0. 0000000000000279130 
9.6251250 0. 0000000000000279135 
9. 6251250 0. 0000000000000279140 
9. 6251250 0. 0000000000000279145 
9.6251250 0. 0000000000000279150 
9.6251250 0. 0000000000000279155 
9. 6251250 0. 0000000000000279160 
= 1.25 mm and «> = 9.8 

FREQUENCY (GHZ) CAPACITANCE (FARAAD) 
9. 7803600 0. 0000000000000318055 
9. 7803600 0. 0000000000000318058 
9. 7803600 0. 0000000000000318061 
9. 7803600 0. 0000000000000318064 
9. 7803600 0. 0000000000000318067 
9. 7803600 0. 0000000000000318070 
9. 7803600 0. 0000000000000318073 
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IMPEDANCE 


0. 
0. 
0. 
0. 
0. 


000205384374 
000119294322 
000033204290 
000052885720 
000138975711 


IMPEDANCE 


oooo°co°o 


ooooo0o°o 


. 000283333299 
. 000175879548 
- 000068425835 
. 000039027837 
. 000146481471 
. 000253935066 
. 000361388622 


IMPEDANCE 


. 000174399228 
- 000123674994 
. 000072950771 
. 000022226558 
. 000028497644 
. 000079221837 
. 000129946020 











“9 


APPENDIX B. PROGRAMS FOR CALCULATION OF RESONANT 


FREQUENCY 


* VARIABLE DEFINITION: 

* Ce — 

* I = REPRESENT THE REGION 1,2,3 IN RESONATOR 

* ALPA = DISCRETE TRANSFORM VARIABLE FOR DOMINANT MODE (x-direction) 
* BETA = DISCRETE TRANSFORM VARIABLE FOR DOMINANT MODE (z-direction) 
* A = A HALF WIDTH OF RESONATOR DEVICE = 1.6 mm 

* D = HEIGHT OF RESONATOR BELOW SUBSTRATE = 0.66 mm 

* T = THE THICKNESS OF SUBSTRATE = 0.255 mm 

* W = THE WIDTH OF STRIP, CAN BE 0.711 mm, 1.0 mm, or 1.25 mm 
* H = THE HEIGHT OF RESONATOR ABOVE SUBSTRATE = 0.66 mm 

* AL = THE LENGTH OF STRIPLINE , ASSUME = 10 mm 

* AKI = RADIAN FREQUENCY * SQRT OF ( PERMITTIVITY*PERMEABILITY) 
* MO = PERMEABILITY IN VACCUM = PI*4. E-7 

* MR = PERMEABILITY IN A SUBSTRATE ASSUME = 1 

* EO = PERMITTIVITY IN VACCUM = (1/(36*PI))*E-9 = 8. 85E-12 

* ER = PERMITIVITY IN A SUBSTRATE, CAN BE 2.2 OR 9.8 

* FREQ = FREQUENCY OPERATION IN GHZ 

* DE = DENOMINATOR OF ZE 

* DE = NUMERATOR OF ZE 

* DH1,DH2,DH = DENOMINATOR OF ZH 

* NH = NUMERATOR OF ZH 

* ZE = IMPEDANCE OF ELECTRIC FIELD, eqn (1g) 

* ZH = IMPEDANCE OF MAGNETIC FIELD, eqn (ih) 

* 22Z,Z22X,2XZ,AND ZXX ARE THE GREEN'S IMPEDANCE FUNCTION (ic,d,e,f) 
* COT1 = CT1, EQN (1i) 

* COT2 = CT2, EQN (1)}) 

* COT3 = CT3, EQN (1k) 

* J1 = EQN (8a) 

* J2 = EQN (8b) 

* JJ1 = Jz1 EQN (6a) 

* JJ4 = Jxl EQN (6b) 

* JZ1 = BASIS FUNCTION OF CURRENT COMPONENT IN Z - DIRECTION 

* JX1 = BASIS FUNCTION OF CURRENT COMPONENT IN X - DIRECTION 

* K11,K12,K21,K22 = MATRIX ELEMENT 

* J1, J3 = ASSUME CURRENT DISTRIBUTION IN X - DIRECTION 

* J2, J& = ASSUME CURRENT DISTRIBUTION IN Z - DIRECTION 

* ALL COMPUTATION IN STANDARD "MKS" 
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FIRST ORDER APPROXIMATION 


THIS IS AN EXAMPLE PROGRAM TO COMPUTE THE VALUE 
OF DETERMINANT OF GREEN FUNCTION IN FIRST ORDER 
APPROXIMATION WITH W = 0.711 mm and ER = 2.2 





MAIN PROGRAM : 


FUNCTION : F1(ALPA,BETA,FREQ) = Z2Z * JZ1 **2 
COMPLEX KK,K11,F1 

INTEGER N,M,L,I 

REAL A,W,D,T,H,MO,E0O,ER,PI,MAGKK,AL,FREQ,AI.PA,BETA 
PARAMETER (PI = 3.141592654, EO = 8. 85E-12) 
PARAMETER (H = 0. 66E-03, T = 0. 255E-03, D = 0. 66E-03) 
PARAMETER (MO = PI*4.E-7, AL = 0.01, A = 1. 6E-03) 
PARAMETER (ER = 2.2 , W = 0.711E-03) 

NOTE : PARAMETER W CAN BE 0.711, 1.0 OR 1.25 mm 

NOTE : PARAMETER ER CAN BE 2.2 OR 9.8 


DO 15 FREQ = 12. 9075,12.9125,3. 34E-4 
KK = CMPLX ( 0.0,0.0 ) 


START SUMMATION IN ALPA FROM N=~30 TO N=+30 


ALPA = (N-0.5)*PI/A 
K11 = CMPLX ( 0.0,0.0 ) 


START SUMMATION IN BETA FROM M=~-3000 TO M=+3000 
DO 35 I = 1,6001 

M= 1 - 3001 

BETA = (M-0.5)*PI/(5. *AL) 

K11 = K11 + F1(ALPA, BETA, FREQ) 

CONTINUE 

KK = KK + K11 

CONTINUE 

MAGKK = CABS((KK)) 


WRITE( 6, 100) FREQ, MAGKK 
FORMAT(” ',3X,F25. 12,3X,F35. 10) 
CONTINUE 

STOP 

END 


FUNCTION SUBPROGRAM F1 


COMPLEX FUNCTION F1(ALPA, BETA, FREQ) 

COMPLEX ZE,ZH,DE,NE,NH,COT1,COT2 ,COT3 ,GM1,GM2 ,GM3 

COMPLEX YE,YH,DH1,DH2,DH,2221,22Z2,22Z2,J21 

REAL G1,G2,G3,AK1,AK2 ,AKK1,AKK2 ,GMM1,GMM2 ,GMM3 ,COTT1,COTT2 ,COTT3 
REAL YYE,YYH,ALPA,BETA,FREQ,AL,H,D,T,W,A,MO,PI,ER,EO 

REAL X,Y,X1,J11,J12,J13,J14,J15,J1,J21,J22,J2,JJ1 
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PARAMETER (PI = 3. 141592654, EO = 8. 85E~-12) 
PARAMETER (H = 0. 66E-03, T = 0. 255E-03, D = 0. 66E-03) 
PARAMETER (MO = PI*4.E-7, AL = 0.01, A = 1.6E-03) 
PARAMETER (W = 0.711E-03, ER = 2.2 ) 

NOTE : PARAMETER W CAN BE 0.711, 1.0 OR 1.25 mm 

NOTE : PARAMETER ER CAN BE 2.2 OR 9.8 


YYE = -1. /(2. E+9*PI*FREQ*E0) 

YE = CMPLX ( 0.0,YYE ) 

YYH = 2. E+9*PI*MO*FREQ 

YH = CMPLX(0. 0,YYH) 

AK1 = MO*EO*( (2. E+9*PI*FREQ)**2) 
AK2 = MO*EO*ER*( (2. E+9*PI*FREQ)**2) 
AKK1 = ALPA%*2 + BETA**2 - AK1 
AKK2 = ALPA**2 + BETA**2 ~ AK2 
G1 = SQRT(ABS(AKK1)) 

G3 = Gl 

IF (AKK1 .LT. 0.0) THEN 

GM1 = CMPLX(0. 0,G1) 


GM3 = GM1 

GMM1 = G1*H 

GMM3 = G3*D 

COTT1 = -1. /TAN(GMM1) 
COTT3 = -1. /TAN(GMM3) 
COT1 = CMPLX(0. 0,COTT1) 
COT3 = CMPLX(0. 0,COTT3) 
ELSE 

GM1 = CMPLX(G1,0. 0) 

GM3 = GM1 

GMM1 = G1*H 

GMM3 = G3*D 

coTT1 1. /TANH(GMM1) 


COTT3 = 1. /TANH(GMM3) 

COT1 = CMPLX(COTT1,0. 0) 
COT3 = CMPLX(COTT3,0. 0) 
END IF 

G2 = SQRT(ABS( AKK2) ) 

IF (AKK2 .LT. 0.0) THEN 
GM2 = CMPLX ( 0.0,G2 ) 
GMM2 = G2*T 

COTT2 = -1. /TAN(GMM2) 

COT2 = CMPLX ( 0.0,COTT2 ) 


ELSE 
GM2 = CMPLX(G2,0. 0) 
GMM2 = G2*T 


COTT2 = 1. /TANH(GMM2) 

COT2 = CMPLX ( COTT2,0.0 ) 

END IF 

DE = COT2*COT3 + COT1*COT3*(GM2/GM1)/ER + COT1*COT2 + (GM3/GM2)*ER 
NE = ((GM2/ER)*COT3 + GM3*COT2)*YE 


NH = GM2*COT2 + GM3*COT3 

DH1 = GM1*GM2*COT1*COT2 + GM1*GM3*COT1*COT3 
DH2 = GM2*GM3*COT2*COT3 + GM2**2 

DH = DH1 + DH2 

ZH = ( NH / DH )*YH 

ZE = NE / DE 
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X1 = -1. /(ALPA**2 + BETA**2) 
ZZZ1 = CMPLX ( X1,0.0 ) 


2222 = (BETA**2)*ZE + (ALPA**2)*ZH 
ZZ2Z = Z22Z1*ZZZ2 

X = ALPA*W/2. 

J11 = 2. *SIN(X)/X 

J12 = 3. /CX**2) 

J13 = COS(X) 

J14 = Jill 


J15 (2. *(1. -COS(X)))/(X**2) 

JL = J11 + J12*(J13 - Ji4 + J15) 
Y = BETA*AL / 2. 

J21 = PI * COS(Y) 

J22 = (PI/2. )¥*2 - Ye*2 

J2 = J21 / J22 

JJ1 = J1*J2 

JZ1 = CMPLX(JJ1,0. 0) 

Fl = 222*JZ1*JZ1 

RETURN 
END 





++ tt 


35 


25 


SECOND ORDER APPROXIMATION 


THIS IS AN EXAMPLE PROGRAM TO COMPUTE THE VALUE 
OF DETERMINANT OF GREEN FUNCTION IN SECOND ORDER 
APPROXIMATION WITH W = 0.711 mm and ER = 2.2 





MAIN PROGRAM : 


F1(ALPA,BETA,FREQ) = ZZZ * JZ1 ** 2 
F2(ALPA,BETA,FREQ) = 2Z2X * JZ1 * JX1 
F3(ALPA,BETA,FREQ) = ZXZ * JZ1 * JX1 
F4(ALPA, BETA,FREQ) = ZXX * JX1 ** 2 


COMPLEX K11,K22,K33,K44,KK1,KK2,KK3,KK4 

COMPLEX KKK1,KKK2 ,KK 

COMPLEX F1,F2,F3,F4 

INTEGER M,N,I,L 

REAL MAGKK,A,W,D,T,H,MO,EO,ER,PI,AL,ALPA, BETA, FREQ 
PARAMETER ( PI=3. 141592654, E0=8. 85E- 12 ) 
PARAMETER ( H=0. 66E-03, T=0. 255E-03, D=0.66E-03 ) 
PARAMETER ( MO=PI*4. E- a AL=0.01, A =1.6E-03 ) 
PARAMETER ( ER = 2.2 , W = 0.711E-03 ) 

NOTE : PAREMETER W CAN BE 0.711, 1.0, OR 1.25 mm 
NOTE : PAREMETER ER CAN BE 2.2 OR 9.8 


DO 15 FREQ = 12. 8582191,12. 8582687,2. 48E-6 
KK1 = CMPLX( 0.0,0.0 ) 
KK2 = CMPLX( 0.0,0.0 ) 
KK3 = CMPLX( 0.0,0.0 ) 
KK4 = CMPLX( 0.0,0.0 ) 
START SUMMATION IN ALPA FROM N=-30 TO N=+30 


DO 25 L= 1, 61 


N=1L- 31 

ALPA = (N-0.5)*PI/A 
K1l = CMPLX ( 0.0,0.0 ) 
K22 = CMPLX ( 0.0,0.0 ) 
K33 = CMPLX ( 0.0,0.0 ) 
K44 = CMPLX ( 0.0,0.0 ) 


DO 35 I= 1, 6001 
M= 1 - 3001 
BETA = (M-0.5)*PI/(5. *AL) 


K1ll = K11 + F1CALPA, BETA, FREQ) 
K22 = K22 + F2(ALPA, BETA, FREQ) 
K33 = K33 + F3(ALPA,BETA,FREQ) 
K44 = K44 + F4CALPA, BETA, FREQ) 
CONTINUE 

KK1 = KK1 + K11 

KK2 = KK2 + K22 

KK3 = KK3 + K33 

KK4 = KK4 + K4&& 

CONTINUE 
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KKK1 = KK1*KK4 
KKK2 = KK2*KK3 
KK = KKK1 - KKK2 
MAGKK = CABS((KK)) 


WRITE(6, 100) FREQ, MAGKK 
100 FORMAT(” ',2X,F16.12,3X,F20. 10) 
15 CONTINUE 

STOP 

END 


* FUNCTIONS SUBPROGRAM F1, F2, F3, F4 


COMPLEX FUNCTION F1(ALPA, BETA, FREQ) 
COMPLEX ZE,ZH,DE,NE,NH,COT1,COT2,COT3 
COMPLEX GM1,GM2,GM3,YE,YH,DH1,DH2,DH 
COMPLEX 2221,2222,22Z,JZ1 
REAL G1,G2,G3,AK1,AK2,AKK1 , AKK2 ,GMM1 ,GMM2 ,GMM3 ,COTT1,COTT2,COTT3 
REAL YYE,YYH,ALPA, BETA,FREQ,AL,H,T,D,W,A,MO,PI,ER,E0 
REAL X,Y,X1,J11,J12,J13,J14,J15,J1,J21,J22,J2,JJ1 
PARAMETER ( PI=3. 141592654, ER=2.2, EO=8.85E-12 ) 
PARAMETER ( H=0. 66E-03, T=0.255E-03, D=0.66E-03 ) 
PARAMETER ( MO=PI*4.E-7, AL=0.01 ) 
PARAMETER ( A=1.6E-03, W=0.711E-03 ) 

* NOTE : PAREMETER W CAN BE 0.711, 1.0, OR 1.25 mm 

* NOTE : PAREMETER ER CAN BE 2.2 OR 9.8 


YYE = -1./ (2. E+9*PI*FREQ*E0) 

YE = CMPLX (0.0,YYE) 

YYH = 2. E+9*PI*MO*FREQ 

YH = CMPLX (0.0,YYH) 

AK1 = MO*EO*( (2. E+9*PI*FREQ)**2) 
AK2 = MO*EO*ER*( (2. E+9*PI*FREQ)**2) 
AKK1 = ALPA**2 + BETA**2 - AK1 
AKK2 = ALPA**2 + BETA**2 - AK2 
G1 = SQRT(ABS(AKK1) ) 

G3 = Gl 

IF (AKK1 .LT. 0.0) THEN 

GM1 = CMPLX ( 0.0,G1 ) 


GM3 = GM1 

GMM1 = G1*H 

GMM3 = G3*D 

CcOTT1 -1. /TAN(GMM1) 


COTT3 = -1. /TAN(GMM3) 


COT1 = CMPLX(0. 0,COTT1) 
COT3 = CMPLX(0. 0,COTT3) 
ELSE 

GM1 = CMPLX ( G1,0.0 ) 

GM3 = GM1 

GMM1 = G1*H 

GMM3 = G3*D 

COTT1 = 1./ TANH(GMM1) 
COTT3 = 1./ TANH(GMM3) 
COT1 = CMPLX ( COTT1,0.0 ) 
COT3 = CMPLX ( COTT3,0.0 ) 
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END IF 

G2=SQRT( ABS( AKK2) ) 

IF (AKK2 .LT. 0.0) THEN 
GM2 = CMPLX ( 0.0,G2 ) 
GMM2 = G2*T 

COTT2 = -1. /TAN(CMM2) 

COT2 = CMPLX ( 0.0,COTT2 ) 


ELSE 
GM2 = CMPLX ( G2,0.0 ) 
GMM2 = G2*T 


COTT2 = 1./ TANH(GMM2) 

COT2 = CMPLX ( COTT2,0.0 ) 

END IF 

DE = COT2*COT3+COT1*COT3*( GM2/GM1) /ER+COT1*COT2+(GM3/GM2)*ER 
NE = ((GM2/ER)*COT3+GM3*COT2 )*YE 

NH= GM2*COT2 + GM3*COT3 

DH1 = GM1*GM2*COT1*COT2 + GM1*GM3*COT1*COT3 

DH2 = GM2*GM3*COT2*COT3 + GM2**2 


DH = DH1+DH2 

ZH = (NH / DH) * YH 

ZE = NE / DE 

X1 = -1./( ALPA**2 + BETA**2 ) 
Z2Z1 = CMPLX ( X1,0.0 ) 

ZZ22 = (BETA**2)*ZE + (ALPA**2)*ZH 


Z2Z = Z2221*2222 
X = ALPA*W/2. 
Y = BETA*AL / 2. 


J11 = 2. *SIN(X)/(X) 

J12 = 3. /(X**2) 

J1i3 = COS(X) 

J14 = Jil 

Ji5 = (2.¥*(1. -COS(X))) / (X**2) 
Jl = J114+J12*(J13-J14+J15) 

J21 = PI*COS(Y) 

J22 = (PI/2. )**2 - (Y¥**2) 

J2 = J21/J22 

JJ1 = J1*J2 


JZ1 = CMPLX ( JJ1,0.0 ) 
Fl = 222*JZ1*JZ1 
RETURN 

END 


COMPLEX FUNCTION F2(ALPA, BETA, FREQ) 

COMPLEX ZE,ZH,DE,NE,NH,COT1,COT2 ,COT3 

COMPLEX GM1,GM2,GM3,YE,YH,DH1,DH2,DH 

COMPLEX 2ZX,JZ1,JX1 

REAL G1,G2,G63,AK1,AK2,AKK1 , AKK2 ,GMM1 ,GMM2 ,GMM3 ,COTT1,COTT2 ,COTT3 
REAL YYE, YYH,ALPA, BETA, FREQ,AL,H,T,D,W,A,MO,PI,ER,EO 
REAL X,Y,J11,J12,J13,J14,J15,J1,J21,J22,J2,JJ1 

REAL J41,J42,J4,J31,J32,J3,JJ4,Z2ZX1 

PARAMETER ( PI=3. 141592654, ER=2.2, EO=8. 85E-12 ) 
PARAMETER ( H=0. 66E-03, T=0. 255E-03, D=0.66E-03 ) 
PARAMETER ( MO=PI*4.E-7, AL=0.01 ) 

PARAMETER ( A=0.16E-02, W=0.711E-03 ) 

NOTE : PAREMETER W CAN BE 0.711, 1.0, OR 1.25 mm 
NOTE : PAREMETER ER CAN BE 2.2 OR 9.8 
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YYE = -1. /(2. E+9*PI*FREQ*E0) 

YE = CMPLX (0.0,YYE) aT) 
YYH = 2. E+9*PI*MO*FREQ 

YH = CMPLX (0.0, YYH) 

AK1 = MO*E0*( (2. E+9*PI*FREQ)**2) 
AK2 = MO*EO*ER*( (2. E+9*PI*FREQ)**2) 
AKK1 = ALPA**2 + BETA**2 - AKi 

AKK2 = ALPA**2 + BETA**2 - AK2 

G1 = SQRT( ABS(AKK1)) 

G3 = Gl 

IF (AKK1 .LT. 0.0) THEN 

GM1 = CMPLX ( 0.0,G1 ) 


GM3 = GM1 
GMM1 = G1*H 
GMM3 = G3*D 


COTT1 = -1. /TAN(GMM1) 
COTT3 = -1. /TAN(GMM3) 
COT1 = CMPLX(0.0,COTT1) 


COT3 = CMPLX(0. 0,COTT3) 
ELSE 

GM1 = CMPLX ( G1,0.0 ) 
GM3 = GM1 

GMM1 = G1*H 

GMM3 = G3*D 


COTT1 = 1. /TANH(GMM1) 
COTT3 = 1. /TANH(GMM3) 

COT1 = CMPLX ( COTT1,0.0 ) 
COT3 = CMPLX ( COTT3,0.0 ) 
END IF 

G2=SQRT( ABS( AKK2) ) 

IF (AKK2 .LT. 0.0) THEN 
GM2 = CMPLX ( 0.0,G2 ) 
GMM2 = G2*T 

COTT2 = -1. /TAN(GMM2) 

COT2 = CMPLX ( 0.0,COTT2 ) 


ELSE 
GM2 = CMPLX ( G2,0.0 ) 
GMM2 = G2*T 


COTT2 = 1. /TANH(GMM2) 
COT2 = CMPLX ( COTT2,0.0 ) 


END IF 

DE = COT2*COT3+COT1*COT3*(GM2/GM1) /ER+COT1*COT2+( GM3/GM2)*ER 
NE = ((GM2/ER)*COT3+GM3*COT2)*YE 

NH = GM2*COT2 + GM3*COT3 

DH1 = GM1*GM2*COT1*COT2 + GM1*GM3*COT1*COT3 

DH2 = GM2*GM3*COT2*COT3 + GM2**2 

DH = DH1+DH2 

ZH = (NH/DH)*YH 

ZE = NE/DE 


ZZX1 = -1. *(ALPA*BETA) /( ALPA**2+BETA** 2) 
ZZX = Z2Z2X1*( ZE-ZH) 

X = ALPA*W/2. 

Y = BETA*AL/2. 

J31 = 2. *PI*SIN(X) 

J32 = X%*2 = PI? 
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J3 = J31/J32 


J41 = COS(Y)/(Y) 

J42 = SINCY)/(Y**2) 

J& = J41-J42 

JI4 = J3*IJ4 

JX1 = CMPLX ( JJ4,0.0 ) 
J11 = 2. *SIN(X)/(X) 

J12 = 3. /(X**2) 

J13 = COS(X) 

J14 = 2. *SIN(X)/(X) 

J15 = (2. *(1. -COSC(X)))/CX**2) 
J1 = J114+J12*(J13-J14+J15) 
J21 = PI*COS(Y) 

J22 = (P1I/2)**2-Y**2 

J2 = J21/J22 

JJ1 = J1*J2 

JZ1 = CMPLX (JJ1,0.0 ) 

F2 = Z2X*JX1*JZ1 

RETURN 

END 


COMPLEX FUNCTION F3(ALPA, BETA, FREQ) 

COMPLEX ZE,ZH,DE,NE,NH,COT1,COT2,COT3 

COMPLEX GM1,GM2,GM3,YE,YH,DH1,DH2,DH 

COMPLEX ZXZ,JZ1,JX1 

REAL G1,G2,G3,AK1,AK2,AKK1,AKK2 ,GMM1,GMM2 ,GMM3,COTT1,COTT2 ,COTT3 
REAL YYE,YYH,ALPA, BETA, FREQ,AL,H,T,D,W,A,MO,PI,ER,E0 
REAL X,Y,J11,J12,J13,J14,J15,J1,J21,J22,J2,JJ1 

REAL J41,J42,J4,J31,J32,J3,JJ4,ZX21 

PARAMETER ( PI=3. 141592654, ER=2.2, E0=8.85E-12 ) 
PARAMETER ( H=0. 66E-03, T=0. 255E-03, D=0.66E-03 ) 
PARAMETER ( MO=PI*4.E-7, AL=0.01 ) 

PARAMETER ( A=1.6E-03, W=0.711E-03 ) 

NOTE : PAREMETER W CAN BE 0.711, 1.0, OR 1.25 mm 
NOTE : PAREMETER ER CAN BE 2.2 OR 9.8 


YYE = -1. /(2. E+9*PI*FREQ*EO) 

YE = CMPLX (0.0,YYE) 

YYH = 2. E+9*PI*MO*FREQ 

YH = CMPLX (0.0,YYH) 

AK1 = MO*EO*( (2. E+9*PI*FREQ)**2) 
AK2 = MO*EO*ER*( (2. E+9*PI*FREQ)**2) 
AKK1 = ALPA**2 + BETA**2 - AK1 
AKK2 = ALPA**2 + BETA**2 - AK2 
G1 = SQRT(ABS( AKK1)) 

G3 = Gl 

IF (AKK1 .LT. 0.0) THEN 

GM1 = CMPLX ( 0.0,G1 ) 


GM3 = GM1 
GMM1 = G1*H 
GMM3 = G3*D 


COTT1 = -1. /TAN(GMM1) 


COTT3 -1. /TANCGMM3) 
COT1 = CMPLX(0.0,COTT1) 
COT3 = CMPLX(0. 0,COTT3) 
ELSE 


GM1 
GM3 


GMM1 
GMM3 
CcOTT1 


CMPLX ( Gi,0.0 ) 
GM1 

= G1*H 

= G3*D 


1. /TANH(GMM1) 


COTT3 = 1. /TANH(GMM3) 


COT1 = CMPLX ( COTT1,0.0 ) 
COT3 = CMPLX ( COTT3,0.0 ) 
ENDIF 


G2=SQRT( ABS( AKK2) ) 
IF (AKK2 .LT. 0.0) THEN 


GM2 


= CMPLX ( 0.0,G2 ) 


GMM2 = G2*T 
COTT2 = -1. /TAN(GMM2) 
COT2 = CMPLX ( 0.0,COTT2 ) 


ELSE 


GM2 


= CMPLX ( G2,0.0 ) 


GMM2 = G2*T 
COTT2 = 1. /TANH(GMM2) 
COT2 = CMPLX ( COTT2,0.0 ) 


END IF 

DE = COT2*COT3+COT1*COT3*(GM2/GM1) /ER+COT1*COT2+( GM3/GM2)*ER 
NE = ((GM2/ER)*COT3+GM3*COT2)*YE 
NH = GM2*COT2 + GM3*COT3 

DH1 = GM1*GM2*COT1*COT2 + GM1*GM3*COT1*COT3 
DH2 = GM2*GM3*COT2*COT3 + GM2%*2 
DH = DH1+DH2 

ZH = (NH / DH) * YH 

ZE = NE / DE 

ZXZ1 = ~1.*(ALPA*BETA) /(ALPA**2 + BETA**2) 
2XZ = ZXZ1*(ZE-ZH) 

X = ALPA*W/2. 

Y = BETA*AL / 2. 

J11 = 2. *SIN(X)/(X) 

J12 = 3. /(X**2) 

J13 = COS(X) 

J14 = Jll 

J15 = (2.0*(1.0-COS(X))) / (X#*2) 
J1 = J114J12*(J13-J3144J15) 

J21 = PI*COS(Y) 

J22 = (PI/2. )**2 - (Y#*2) 

J2 = J21/J22 

JJ1 = Ji*J2 

JZ1 = CMPLX ( JJ1,0.0 ) 

J31 = 2.*PI*SIN(X) 

J32 = Xwe2 = PI**?2 

J3 = J31 / 332 

J41 = COS(Y)/(Y) 

J42 = SINC(Y)/(Y#*2) 

J4 = J41-J42 

JI4 = I3*I4 

JX1 = CMPLX ( JJ4,0.0 ) 


F3 = ZXZ*JZi*JX1 
RETURN 


END 
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COMPLEX FUNCTION F4(ALPA, BETA, FREQ) 
COMPLEX ZE,ZH,DE,NE,NH,COT1,COT2,COT3 
COMPLEX GM1,GM2,GM3,YE,YH,DH1,DH2,DH 
COMPLEX ZXX1,ZXX2,ZXX,JX1 


REAL G1,G2,G3,AK1,AK2,AKK1,AKK2 ,GMM1,GMM2 ,GMM3,COTT1,COTT2,COTT3 
REAL YYE,YYH,ALPA,BETA,FREQ,AL,H,T,D,W,A,MO,PI,ER,EO 


REAL X,Y,X1,J31,J32,J3,J41,J42,J4,JJ54 


PARAMETER ( PI=3. 141592654, ER=2.2, E0=8.85E~-12 ) 
PARAMETER ( H=0. 66E-03, T=0.255E-03, D=0.66E-03 ) 


PARAMETER ( MO=PI*4.E-7, AL=0.01 ) 
PARAMETER ( A=1.6E-03, W=0.711E-03 ) 


NOTE : PAREMETER W CAN BE 0.711, 1.0, OR 1.25 mm 


NOTE : PAREMETER ER CAN BE 2.2 OR 9.8 


YYE = -1. /(2. E+9*PI*FREQ*E0) 

YE = CMPLX (0. 0,YYE) 

YYH = 2. E+9*PI*MO*FREQ 

YH = CMPLX (0.0, YYH) 

AK1 = MO*EO*( (2. E+9*PI*FREQ)**2) 
AK2 = MO*EO*ER*( (2. E+9*PI*FREQ)**2) 
AKK1 = ALPA**2 + DETA**2 - AKI 
AKK2 = ALPA**2 + BETA**2 - AK2 
G1 = SQRT(ABS(AKK1)) 

G3 = Gl 

IF (AKK1 .LT. 0.0) THEN 

GM1 = CMPLX ( 0.0,G1 ) 


GM3 = GM1 

GMM1 = G1*H 

GMM3 = G3*D 

COTT1 = -1. /TAN(GMM1) 


COTT3 = -1. /TAN(GMM3) 


COT1 = CMPLX(0. 0,COTT1) 
COT3 = CMPLX(0. 0,COTT3) 
ELSE 

GM1 = CMPLX ( G1,0.0 ) 
GM3 = GM1 

GMM1 = G1*H 

GMM3 = G3*D 


COTT1 = 1. /TANH(GMM1) 
= 1. /TANH(GMM3) 
CMPLX ( COTT1,0. 


cOT1 0) 
CMPLX ( COTT3,0.0 ) 


coT3 
ENDIF 
G2=SQRT( ABS( AKK2) ) 

IF (AKK2 .LT. 0.0) THEN 
GM2 = CMPLX ( 0.0,G2 ) 
GMM2 = G2*T 

COTT2 = -1. /TAN(GMM2) 

COT2 = CMPLX ( 0.0,COTT2 ) 


ELSE 
GM2 = CMPLX ( G2,0.0 ) 
GMM2 = G2*T 


COTT2 = 1. /TANH(GMM2) 
COT2 = CMPLX ( COTT2,0.0 ) 
END IF 


DE = COT2*COT3+COT1*COT3*( GM2/GM1) /ER+COT1*COT2+( GM3/GM2)*ER 
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NE = ((GM2/ER)*COT3+GM3*COT2)*YE 
NH = GM2*COT2 + GM3*COT3 
DH1 = GM1*GM2*COT1*COT2 + GM1*GM3*COT1*COT3 


DH2 = GM2*GM3*COT2*COT3 + GM2**2 
DH = DH1+DH2 

ZH = (NH / DH) * YH 

ZE = NE / DE 

Xi = -1. /(ALPA**2 + BETA**2) 
ZXX1 = CMPLX (X1,0.0) 


ZXX2 = (ALPA**2)*ZE + (BETA**2)*ZH 
ZXX = ZXX1*ZXX2 

X = ALPA*W/2. 

Y = BETA*AL/2. 

J31 = 2. *PI*SIN(X) 

J32 = X**2 - PI #2 

J3 = J31 / J32 

J41 = COS(Y)/(Y) 


J42 = SINCY)/CY**2) 
J4 = J41-J42 
JI&4 = J3*I4 


JX1 = CMPLX ( JJ4,0.0 ) 
F4 = ZXX*JX1*JX1 

RETURN 

END 
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APPENDIX C. PROGRAMS FOR CALCULATION OF 


PROPAGATION CONSTANT 


VARIABLE DEFINITION: 


T 
ALPA 
BETA 


> 
co 
Unit od bet tb ton ea 


REPRESENT THE REGION 1,2,3 IN RESONATOR 

DISCRETE TRANSFORM VARIABLE FOR DOMINANT MODE (x-direction) 
DISCRETE TRANSFORM VARIABLE FOR DOMINANT MODE (2-direction) 
A HALF WIDTH OF RESONATOR DEVICE = 1.6 mm 

HEIGHT OF RESONATOR BELOW SUBSTRATE = 0.66 mm 

THE THICKNESS OF SUBSTRATE = 0.255 mm 

THE WIDTH OF STRIP, CAN BE 0.711 mm, 1.0 mm, or 1.25 mm 
THE HEIGHT OF RESONATOR ABOVE SUBSTRATE = 0.66 mm 

THE LENGTH OF STRIPLINE , ASSUME = 10 mm 

RADIAN FREQUENCY * SQRT OF (PERMITTIVITY*PERMEABILITY) 
PERMEABILITY IN VACCUM = PI*4.E-7 

PERMEABILITY IN A SUBSTRATE ASSUME = 

PERMITTIVITY IN VACCUM = (1/(36*PI))*E-9 = 8. 85E~-12 
PERMITIVITY IN A SUBSTRATE, CAN BE 2.2 OR 9.8 

FREQUENCY OPERATION IN GHZ 

DENOMINATOR OF ZE 

NUMERATOR OF ZE 


DH1,DH2,DH = DENOMINATOR OF ZH 


NH = 


N 
3] 
woul 


NUMERATOR OF ZH 
IMPEDANCE OF ELECTRIC FIELD, eqn (1g) 
IMPEDANCE OF MAGNETIC FIELD, eqn (ih) 


22Z,2Z2X,ZXZ,AND ZXX ARE THE GREEN'S TNPEDANCE FUNCTION (1c,d,e,f) 


COT1 
COT2 
COoT3 
Jl 
J2 
JJ1 
JI4 
JZ1 
JX1 


Vine daw da 


CT1, EQN (14) 

CT2, EQN (1)) 

CT3, EQN (1k) 

EQN (8a) 

EQN (8b) 

Jz1 EQN (6a) 

Jx1 EQN (6b) 

BASIS FUNCTION OF CURRENT COMPONENT IN Z - DIRECTION 
= BASIS FUNCTION OF CURRENT COMPONENT IN X - DIRECTION 


Kll, K12, K21,K22 = MATRIX ELEMENT 


J1, J3 
J2, J4 


= ASSUME CURRENT DISTRIBUTION IN X - DIRECTION 
= ASSUME CURRENT DISTRIBUTION IN Z - DIRECTION 


ALL COMPUTATION IN STANDARD "MKS" 
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FIRST ORDER APPROXIMATION 


THIS IS AN EXAMPLE PROGRAM TO COMPUTE CONSTANT PROPAGATION 
IN FIRST ORDER APPROXIMATION WITH W = 0.711 AND ER = 2.2 





MAIN PROGRAM : 


FUNCTION : F1(ALPA,BETA,FREQ) = 222 * JZ1 ** 2 

COMPLEX F1,K11 

REAL BO,B00,MAGKK,ALPA, BETA,MO 

REAL A,D,T,H,PI,ER,W,FREQ,OMG,EO 

INTEGER N,L 

PARAMETER ( PI=3. 141592654, ER=2.2, E0=8.85E-12, H=0.66E-03 ) 
PARAMETER ( T=0. 255E-03, D=0.66E-03, MO=PI*4.E-7, A=0.16E-02 ) 
PARAMETER ( W=0.711E-3 ) 

NOTE : PARAMETER W CAN BE 0.711 MM, 1.0 MM, 1.25 MM 

NOTE : PARAMETER ER CAN BE 2.2 OR 9.8 


FREQ = 12. 9095783 
OMG = FREQ*PI*2, E+9 
BO = (OMG/3. E08)*SQRT( (ER+1)/2) 


START SWEEP THE BETA VALUE FROM 84% TO 90% 
DO 20 BOO = 0. 84,0.9,1. 2E-3 

BETA = BO*BOO 

Kl1 = .0.0,0.0) 


START SUMMATION IN ALPA FROM N=-40 TO N=+40 
DO 10 L=1, 81 

N=L- 41 

ALPA = (N - 0.5)*PI/A 

K1li = Kil + F1(ALPA, BETA, FREQ) 
CONTINUE 

MAGKK = CABS((K11)) 

WRITE(6,100) FREQ, BETA,MAGKK 

FORMAT(' ',3X,F14.9,3X,F18.9,3X,F15.9) 
CONTINUE 

STOP 

END 


FUNCTION SUBPROGRAM F1 


COMPLEX FUNCTION F1(ALPA, BETA, FREQ) 

COMPLEX GM1,GM2,GM3,COT1,COT2 ,COT3,DE,ZE,DH,ZH,NE,NH,DH1,DH2,YE 
COMPLEX JZ1,2221,2222,222, YH 

REAL A,D,T,H,PI,ER,W,FREQ,OMG,E0,G1,G2,G3,COTT1,COTT2,COTT3,YYE 
REAL YYH,GMM1,GMM2,GMM3 ,ALPA,AKK1,AKK2 ,X1,X 

REAL J11,J12,J13,J14,J15,J1,MO,AK1,AK2,BETA 

PARAMETER ( PI=3. 141592654, ER=2.2, EO=8.85E-12, H=0.66E-03 ) 
PARAMETER ( T=0. 255E-03, D=0.66E-03, MO=PI*4.E-7, A=0,16E-02 ) 








PARAMETER ( W=0.711E-3 ) 
NOTE : PARAMETER W CAN BE 0.711 MM, 1.0 MM, 1.25 MM 


NOTE : PARAMETER ER CAN BE 2.2 OR 9.8 fe 
OMG = FREQ*PI*2. E+9 

YYE = -1. /C(OMG*EO) 

YE = CMPLX (0.0,YYE) 

YYH = OMG*MO 

YH = CMPLX (0.0, YYH) 

AK1 = (OMG**2)*MO*EO 


AK2 = (OMG**2)*MO*EO*ER 
AKK1 = ALPA**2 + BETA**2 - AKI 
AKK2 = ALPA**2 + BETA**2 - AK2 
G1 = SQRT(ABS( AKK1)) 

G3 = Gl 

IF (AKK1 .LT. 0.0) THEN 
GM1 = CMPLX ( 0.0,G1 ) 
GM3 = GM1 

GMM1 = G1*H 

GMM3 = G3*D 

CcOTT1 -1. /TAN(GMM1) 
COTT3 -1. /TAN(GMM3) 
COT1 = CMPLX ( 0.0,COTT1 


a 


COT3 = CMPLX ( 0.0,COTT3 
ELSE 

GM1 = CMPLX ( G1,0.0 ) 
GM3 = GM1 

GMM1 = G1*H 

GMM3 = G1*D 

COTT1 = 1. /TANH(GMM1) 
COTT3 = 1. /TANH(GMM3) 


COT1 = CMPLX ( COTT1,0.0 ) 
COT3 = CMPLX ( COTT3,0.0 ) 
END IF 

G2 = SQRT( ABS( AKK2) ) 

IF (AKK2 .LT. 0.0) THEN 
GM2 = CMPLX ( 0.0,G2 ) 
GMM2 = G2*T 

COTT2 = -1. /TAN(GMM2) 

COT2 = CMPLX ( 0.0,COTT2 ) 


ELSE 
GM2 = CMPLX ( G2,0.0 ) 
GMM2 = G2*T 


COTT2 = 1. /TANH(GMM2) 

COT2 = CMPLX ( COTT2,0.0 ) 

END IF 

DE = COT2*COT3+COT1*COT3*( GM2/GM1) /ER+COT1*COT2+( GM3/GM2)*ER 


NE = (GM2*COT3/ER+GM3*COT2)*YE 

ZE = NE/DE 

NH = GM2*COT2 + GM3*COT3 

DH1 = GM1*GM2*COT1*COT2 + GM1*GM3*COT1*COT3 
DH2 = GM2*GM3*COT2*COT3 + GM2**2 

DH = DH1 + DH2 

ZH = (NH/DH)*YH 

Xl = -1. / (ALPA**2 + BETA**2) 


Z2Z1 = CMPLX (X1,0.0) 








Z2Z2 = (BETA**2)*ZE + (ALPA**2)*ZH 
Z2Z = Z22Z1*Z222 


X = ALPA*W/2. 

J11 = 2. *SIN(X)/(X) 
J12 = 3. /(X**2) 

J13 = COS(X) 

Ji4 = Jl 


J15 = 2.%*(1. 0-COS(X))/(X**2) 
Jl = J11+J12*(J13-J14+J15) 
JZ1 = CMPLX ( J1,0.0 ) 

Fl = Z22*JZ1**2 

RETURN 

END 
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SECOND ORDER APPROXIMATION 


THIS IS AN EXAMPLE PROGRAM TO COMPUTE CONSTANT PROPAGATION 
IN SECOND ORDER APPROXIMATION WITH W = 0.711 AND ER = 2.2 


MAIN PROGRAM : 


F1(ALPA, BETA, FREQ) ZZZ * JZ1 ** 2 


F2(ALPA, BETA,FREQ) = ZZX * JZ1 * JX1 
F3(ALPA,BETA,FREQ) = ZXZ * JZ1 * JX1 
F4(ALPA,BETA,FREQ) = ZXX * JX1 ** 2 


COMPLEX F1,F2,F3,F4 

COMPLEX KK,K11,K22,K33,K44 

COMPLEX KKK1,KKK2 

REAL MAGKK,B00,B0,ALPA, BETA,MO 

REAL A,D,T,H,PI,ER,W,FREQ,OMG,EO 

INTEGER N,L 

PARAMETER ( PI=3. 141592654, EO=8.85E-12, H=0.66E-03 ) 
PARAMETER ( T=0. 255E-03, D=0.66E-03, MO=PI*4.E~7, A=0.16E-02 ) 
PARAMETER ( W=0.711 ,ER=2.2 ) 

NOTE : PARAMETER W CAN BE 0.711 , 1.0 , OR 1.25 
NOTE : PARAMETER ER CAN BE 2.2 OR 9.8 


FREQ = 12. 858241 
OMG = FREQ*PI*2. E+9 
BO = (OMG/3. EG8)*SQRT( (ER+1}/2) 


START SWEEP THE BETA VALUE FROM 88.2% TO 89. 7% 
DO 20 BOO = O. 882,0. 897,7. 5E-6 
BETA = BO*BOO 


K1l = (0.0,0. 0) 
K22 = (0.0,0.0) 
K33 = (0.0,0.0) 
K44 = (0.0,0. 0) 


START SUMMATION IN ALPA FROM N=-40 TO N=+40 
DO 10 L=1, 81 


ALPA = (N - 0.5)*PI/A 


Kil = Kill + FICALPA, BETA, FREQ) 
K22 = K22 + F2(ALPA, BETA, FREQ) 
K33 = K33 + F3(ALPA, BETA, FREQ) 
K44 = K44 + F4(ALPA, BETA, FREQ) 
CONTINUE 


KKK1 = K11*K44 

KKK2 = K22*K33 

KK = KKK1-KKK2 

MAGKK = CABS((KK)) 

WRITE(6,100) FREQ, BETA,MAGKK 

FORMAT(' ',3X,F14.9,3X,F18.9,3X,F15. 9) 
CONTINUE 
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STOP 
END 


FUNCTIONS SUBPROGRAM F1, F2, F3, F4 








COMPLEX FUNCTION F1(ALPA, BETA, FREQ) 

COMPLEX GM1,GM2,GM3,COT1,COT2,COT3,DE,ZE,DH,ZH,NE,NH,DH1,DH2,YE 
COMPLEX YH,JZ1,2221,2222,222 

REAL A,D,T,H,PI,ER,W,FREQ,OMG,E0,G1,G62,G3,COTT1,COTT2 ,COTT3,YYE 
REAL YYH,GMM1,GMM2 ,GMM3,ALPA,AKK1,AKK2 ,X1,X 

REAL J11,J12,J13,J14,J15,J1,M0O,AK1,AK2, BETA 

PARAMETER ( PI=3. 141592654, E0=8.85E-12, H=0.66E-03 ) 

PARAMETER ( T=0. 255E-03, D=0.66E-03, MO=PI*4.E-7, A=0.16E-02 ) 
PARAMETER ( W=0.711 ,ER=2.2 ) 

NOTE : PARAMETER W CAN BE 0.711 , 1.0 , OR 1.25 

NOTE : PARAMETER ER CAN BE 2.2 OR 9.8 


OMG = FREQ*PI*2. E+9 
YYE = -1. /COMG*EO) 
YE = CMPLX (0.0,YYE) 
YYH = OMG*MO 

YH = CMPLX (0. 0,YYH) 
AK1 = (COMG**2)*MO*E0O 


AK2 = (OMG**2)*MO*EO*ER 

AKK1 = ALPA**2 + BETA**2 - AK1 
AKK2 = ALPA**2 + BETA**2 - AK2 
G1 = SQRT(ABS(AKK1)) 

G3 = Gl 

IF (AKK1 .LT. 0.0) THEN 

GM1 = CMPLX ( 0.0,G1 ) 


GM3 = GMl 

GMM1 = G1*H 

GMM3 = G3*D 

COTT1 = -1. /TAN(GMM1) 


COTT3 = -1. /TAN(GMM3) 
COT1 = CMPLX ( 0.0,COTT1 ) 


COT3 = CMPLX ( 0.0,COTT3 ) 
ELSE 

GM1 = CMPLX ( G1,0.0 ) 

GM3 = GM1 

GMM1 = G1*H 

GMM3 = G1*D 

COTT1 = 1. /TANH(GMM1) 
COTT3 = 1. /TANH(GMM3) 


COT1 = CMPLX ( COTT1,0.0 ) 
COT3 = CMPLX ( COTT3,0.0 ) 
END IF 

G2 = SQRT(ABS(AKK2) ) 

IF (AKK2 .LT. 0.0) THEN 
GM2 = CMPLX ( 0.0,G2 ) 
GMM2 = G2*T 

COTT2 = -1. /TAN(GMM2) 

COT2 = CMPLX ( 0.0,COTT2 ) 
ELSE 

GM2 = CMPLX ( G2,0.0 ) 








GMM2 = G2*T 
COTT2 = 1. /TANH(GMM2) 
COT2 = CMPLX ( COTT2,0.0 } 


END IF 
DE = COT2*COT3+COT1*COT3*( GM2/GM1) /ER+COT1*COT2+( GM3/GM2)*ER 
NE = (GM2*COT3/ER+GM3*COT2)*YE 
ZE = NE/DE 
NH = GM2*COT2 + GM3*COT3 
DH1 = GM1*GM2*CCT1i*COT2 + GM1*GM3*COT1*COT3 
DH2 = GM2*GM3*COT2*COT3 + GM2**2 
re = DH1 + DH2 
= (NH/DH)*YH 
rai = -1. / (ALPA**2 + BETA**2) 
2221 = CMPLX (X1,0.0) 
ZZZ2 = (BETA**2)*ZE + (ALPA**2)*ZH 
ZZZ = Z2Z22Z1*ZZZ2 
X = ALPA*W/2. 
J1i = 2. *SIN(X)/(X) 
J12 = 3. /(X**2) 
J13 = COS(X) 
J14 = J11 
J15 = 2.*(1. 0-COS(X))/(X**2) 


= J114+J12*(J13-J14+J15) 
JZ1 = CMPLX ( J1,0.0 ) 
Fl = 222*JZ1**2 
RETURN 
END 


COMPLEX FUNCTION F2(ALPA, BETA, FREQ) 

COMPLEX GM1,GM2,GM3,COT1,COT2,COT3,DE,ZE,DH,ZH,NE,NH,DH1,DH2, YE 
COMPLEX YH,JZ1,JX1,22X 

REAL A,D,T,H,PI,ER,W,FREQ,OMG,E0,G1,G2,G3,COTT1,COTT2,COTT3,YYE 
REAL YYH,GMM1,GMM2,GMM3,ALPA,AKK1,AKK2,X,ZZX1,BETA 

REAL J31,J32,J3,J11,J12,J13,J14,J15,J1,M0,AK1,AK2 

PARAMETER ( PI=3. 141592654, EO=8. 85E-12, H=0.66E-03 ) 

PARAMETER ( T=0. 255E-03, D=0.66E-03, MO=PI*4.E-7, A=0.16E-02 ) 
PARAMETER ( W=0.711 ,ER=2.2 ) 

NOTE : PARAMETER W CAN BE 0.711 , 1.0 , OR 1.25 
NOTE : PARAMETER ER CAN BE 2.2 OR 9.8 

OMG = FREQ*PI*2. E+9 

YYE = -1. /(OMG*EO) 

YE = CMPLX (0.0,YYE) 


YYH = OMG*MO 
oe CMPLX (0.0, YYH) 
(OMG#*2)*MO*EO 


AK? = (OMG**2)*MO*EO*ER 
AKK1 = ALPA**2 + BETA**2 - AK1 
AKK2 = ALPA%*2 + BETA**2 - AK2 
G1 = SQRT(ABS(AKK1) ) 

G3 = Gl 

IF (AKK1 .LT. 0.0) THEN 

GM1 = CMPLX ( 0.0,G1 ) 


GM3 = GM1 
GMM1 = G1*H 
GMM3 = G3*D 








COTT1 = -1. /TAN(GMM1) 
COTT3 = ~1. /TAN(GMM3) 

COT1 = CMPLX ( 0.0,COTT1 ) 
COT3 = CMPLX ( 0.0,COTT3 ) 
ELSE 

GM1 = CMPLX ( G1,0.0 ) 

GM3 = GMi 

GMM1 = G1*H 

GMM3 = G1*D 


COTT1 = 1. /TANH(GMM1) 
COTT3 = 1. /TANH(GMM3) 
COT1 = CMPLX ( COTT1,0.0 ) 
COT3 = CMPLX ( COTT3,0.0 ) 
END IF 

= SQRT( ABS( AKK2) ) 
IF (AKK2 .LT. 0.0) THEN 
GM2 = CMPLX ( 0.0,G2 ) 
GMM2 = G2*T 
COTT2 = -1. /TAN(GMM2) 
COT2 = CMPLX ( 0.0,COTT2 ) 


ELSE 
GM2 = CMPLX ( G2,0.0 ) 
GMM2 = G2*T 


COTT2 = 1. /TANH(GMM2) 
COT2 = CMPLX ( COTT2,0.0 ) 


END IF 

DE = COT2*COT3+COT1*COT3*( GM2/GM1) /ER+COT1*COT2+( GM3/GM2)*ER 
NE = (GM2*COT3/ER+GM3*COT2)*YE 

ZE = NE/DE 

NH = GM2*COT2 + GM3*COT3 


DH1 = GM1*GM2*COT1*COT2 + GM1*GM3*COT1*COT3 
DH2 = GM2*GM3*COT2*COT3 + GM2**2 

oe = DH1 + DH2 

= (NH/DH)*YH 

1 = -1. *(ALPA*BETA) /( ALPA**2+BETA**2 ) 
ZZX = Z2X1*( ZE-ZH) 

X = ALPA*W/2. 


J11 = 2. *SINCX)/(X) 

J12 = 3. /(X**2) 

J13 = COS(X) 

J14 = Jl. 

JIS = 2.*(1. 0-COS(X) ) /(CX**2) 


Jl = J11+J12*(J13-J14+J15) 
J31 = 2.*PI*SIN(X) 
J32 = X**2 - PI*#2 
J3 = J31/J332 
JZ1 = CMPLX ( J1,0.0 ) 
JX1 = CMPLX(J3,0. 0) 
= Z2Z2X*JZ1*IX1 
RETURN 
END 


COMPLEX FUNCTION F3(ALPA, BETA, FREQ) 

COMPLEX GM1,GM2,GM3,COT1,COT2,COT3,DE,ZE,DH,2ZH,NE,NH,DH1,DH2,YE 
COMPLEX YH,JZ1,JX1,ZXZ 

REAL A,D,T,H,PI,ER,W,FREQ,OMG,E0,G1,G62,G63,COTT1,COTT2,COTT3, YYE 
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REAL YYH,GMM1,GMM2,GMM3,ALPA,AKK1,AKK2,X,BETA 

REAL J11,J12,J13,J14,J15,J1,J31,J32,J3,M0,AK1,AK2 

PARAMETER ( PI=3. 141592654, E0=8. 85E-12, H=0.66E-03 ) 
PARAMETER ( T=0. 255E-03, D=0. 66E-03, MO=PI*4.E-7, A=0.16E-02 ) 
PARAMETER ( W=0.711 ,ER=2.2 ) 

NOTE : PARAMETER W CAN BE 0.711 , 1.0 , OR 1.25 

NOTE : PARAMETER ER CAN BE 2.2 OR 9.8 


OMG = FREQ*PI*2. E+9 

YYE = -1. /(OMG*EO) 

YE = CMPLX (0.0,YYE) 
YYH = OMG*MO 

YH = CMPLX (0. 0,YYH) 
AK1 = (OMG**2)*MO*EO 
AK2 = (OMG**2)*MO*EO*ER 


AKK1 = ALPA**2 + BETA**2 - AK1 
AKK2 = ALPA**2 + BETA**2 - AK2 
G1 = SQRT( ABS( AKK1) ) 

G3 = Gl 

IF (AKK1 .LT. 0.0) THEN 

GM1 = CMPLX ( 0.0,G1 ) 


GM3 = GM1 

GMM1 = G1*H 

GMM3 = G3*D 

COTT1 -1. /TAN(GMM1) 


COTT3 = -1. /TAN(GMM3) 

COT1 = CMPLX ( 0.0,COTT1 ) 
COT3 = CMPLX ( 0.0,COTT3 ) 
ELSE 

GM1 = CMPLX ( G1,0.0 ) 

GM3 = GM1 

GMM1 = G1*H 

GMM3 = G1*D 


COTT1 = 1. /TANH(GMM1) 
COTT3 = 1. /TANH(GMM3) 

COT1 = CMPLX ( COTT1,0.0 ) 
COT3 = CMPLX ( COTT3,0.0 ) 
END IF 

G2 = SQRT(ABS( AKK2) ) 

IF (AKK2 .LT. 0.0) THEN 
GM2 = CMPLX ( 0.0,G2 ) 
GMM2 = G2*T 

COTT2 = -1. /TAN(GMM2) 

COT2 = CMPLX ( 0.0,COTT2 ) 


ELSE 
GM2 = CMPLX ( G2,0.0 ) 
GMM2 = G2*T 


COTT2 = 1. /TANH(GMM2) 

COT2 = CMPLX ( COTT2,0.0 ) 

END IF 

DE = COT2*COT3+COT1*COT3*( GM2/GM1) /ER+SOT1*COT2+( GM3/GM2)*ER 


NE = (GM2*COT3/ER+GM3*COT2)*YE 

ZE = NE/DE 

NH = GM2*COT2 + GM3*COT3 

DH1 = GM1*GM2*COT1*COT2 + GM1*GM3*COT1*COT3 


DH2 = GM2*GM3*COT2*COT3 + GM2**2 
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DH = DH1 + DH2 
ZH = (NH/DH)*YH 
Z2XZ = -1.*(ALPA*BETA) /( ALPA**2+BETA**2) ag 


Z2XZ = ZX2Z*( ZE-ZH) 
X = ALPA*W/2. 

J31 = 2. *PI*SIN(X) 
J32 = X**2 - PI ee? 
J3 = J31/J32 


JX1 = CMPLX(J3,0. 0) 
J11 = 2. *SIN(X)/(X) 

J12 = 3. /(X**2) 

J13 = COS(X) 

J14 = Jl 

J15 = 2.*(1. 0-COS(X)) /CX**2) 


J1 = J114+J12*(J13-5144+J15) 
JZ1 = CMPLX ( J1,0.0 ) 

F3 = 2X2*JZ1*JX1 

RETURN 

END 


COMPLEX FUNCTION F4(ALPA, BETA, FREQ) 

COMPLEX GM1,GM2,GM3,COT1,COT2,COT3,DE,2°,DH,ZH,NE,NH,DH1,DH2,YE 
COMPLEX YH,JX1,ZXX1, Z2XX2 , ZXX 

REAL A,D,T,H,PI,ER,W,FREQ,OMG,E0,G1,G2,G3,COTT1,COTT2,COTT3,YYE 
REAL YYH,GMM1,GMM2,GMM3,ALPA,AKK1, AKK2,AK1,AK2 

REAL J31,J32,J3,M0,X,X1,BETA 

PARAMETER ( PI=3. 141592654, EO=8. 85E-12, H=0.66E-03 ) 

PARAMETER ( T=0. 255E-03, D=0.66E-03, MO=PI*4.E-7, A=0.16E-02 ) 
PARAMETER ( W=0.711 ,ER=2.2 ) 

NOTE : PARAMETER W CAN BE 0.711 , 1.0 , OR 1.25 

NOTE : PARAMETER ER CAN BE 2.2 OR 9.8 


FREQ*PI*2. E+9 

-1. /COMG*E0) 

YE = CMPLX (0.0,YYE) 

YYH = OMG*MO 

YH = CMPLX (0.0,YYH) 

AK1 = (OMG**%2)*MO*EO 

AK2 = (OMG**2)*MO*EO*ER 

AKK1 = ALPA**2 + BETA**2 - AK1 
AKK2 = ALPA*®*2 + BETA**2 - AK2 
G1 = SQRT(ABS(AKK1) ) 

G3 = Gl 

IF (AKK1 .LT. 0.0) THEN 

GM1 = CMPLX ( 0.0,G1 )} 


3 
tou 


GM3 = GM1 

GMM1 = G1*H 

GMM3 = G3*D 

coTT1 ~1. /TAN(GMM1) 


COTT3 = -1. /TAN(GMM3) 
COT1 = CMPLX ( 0.0,COTT1 


ed ee 


COT3 = CMPLX ( 0.0,COTT3 
ELSE 

GM1 = CMPLX ( G1,0.0 ) 
GM3 = GM1l 


GMM1 = G1*H 








GMM3 = G1*D 
COTT1 = 1. /TANH(GMM1) 
COTT3 = 1. /TANH(GMM3) 

= CMPLX ( COTT1,0.0 ) 
COT3 = CMPLX ( COTT3,0.0 ) 


G2 = SQRT( ABS( AKK2)) 

IF (AKK2 .LT. 0.0) THEN 
GM2 = CMPLX ( 0.0,G2 ) 
GMM2 = G2*T 

COTT2 = -1. /TAN(GMM2) 

COT2 = CMPLX ( 0.0,COTT2 ) 


ELSE 
GM2 = CMPLX ( G2,0.0 ) 
GMM2 = G2*T 


COTT2 = 1. /TANH(GMM2) 
COT2 = CMPLX ( COTT2,0.0 ) 


END IF 

DE = COT2*COT3+COT1*COT3*(GM2/GM1) /ER+COT1*COT2+( GM3/GM2)*ER 
NE = (GM2*COT3/ER+GM3*COT2)*YE 

ZE = NE/DE 

NH = GM2*COT2 + GM3*COT3 


DH1 = GM1*GM2*COT1*COT2 + GM1*GM3*COT1*COT3 
= GM2*GM3*COT2*COT3 + GM2**2 
= DH1 + DH2 
ZH = (NH/DH)*YH 
X1 = -1./(ALPA**2 + BETA**2) 
= CMPLX(X1,0. 0) 
ZXX2 = (ALPA**2)*ZE + (BETA**2)*ZH 
ZXX = ZXX1*ZXX2 
= ALPA*W/2. 
J31 = 2. *PI*SIN(X) 
J32 = X¥*#2 - Pw? 
J3 = J31/J32 
JX1 = CMPLX(J3,0. 0) 
F4 = ZXX*JX1**2 
RETURN 
END 














APPENDIX D. PROGRAMS FOR CALCULATION OF FRINGING 


A. 


t+ + t 


100 


+t + 


100 


CAPACITANCE 


TRANSMISSION LINE CIRCUIT MODEL 


THIS IS AN EXAMPLE PROGRAM TO COMPUTE THE FRINGING CAPACITANCE 
IN SECOND ORDER APPROXIMATION WITH W = 0.711 MM AND ER = 2.2 


COMPLEX ZIN 

REAL MAGZIN,OMG, ZO, BETA,AL,C1,C,FREQ,T,2Z,PI,ZZIN 
PARAMETER ( PI=3. 141592654, Z0=94.088, AL=0.005 ) 
PARAMETER ( BETA = 301.921744 ) 

NOTE : FILL IN PARAMETERS BETA AND ZO FROM TABLE 1 OR TABLE 2. 
DO 10 Cl = 0. 8058,0. 8062,1.E-6 

C = C1*1.E-14 

DO 15 FREQ = 12. 858241,12. 858241 

OMG = 2*PI*FREQ*1.E9 

T = TAN(BETA*AL) 

ZZ = OMG*C*ZO0 

ZZIN = ~-1. *(1-ZZ*T)*Z0/(ZZ+T) 

ZIN = CMPLX(0. 0,2ZIN) 

MAGZIN = CABS((ZIN)) 

WRITE(6 , 100) FREQ,C ,MAGZIN 

FORMAT(' ',2X,F16.12,2X,F25. 23,2X,F14. 12) 
CONTINUE 

CONTINUE 

END 


END EFFECT EXTENSION MODEL 


THIS IS AN EXAMPLE PROGRAM TO COMPUTE THE FRINGING CAPACITANCE 
IN SECOND ORDER APPROXIMATION WITH W = 0.711 MM AND ER = 2.2 


REAL FO,BETA,LL,L,PI,DELTA,CO,OMG,Z0,CP 
PARAMETER ( Z0=94.088, L=0.01, PI=3. 141592654 ) 
PARAMETER ( FO0=12. 858241, BETA=301. 921744 ) 
NOTE : PARAMETERS FO, BETA, AND ZO CAN VARIES DEPEND ON THE CHOICE 
OF APPROXIMATION 

OMG) = 2.*PI*FO*1. E9 

LL = PI/BETA 

DELTA = (LL-L)/2. 

CO = BETA/(OMG*Z0) 

CP = CO*DELTA 

WRITE(6 , 100) FO, BETA,CP 

FORMAT(” ',3X,F10.7,5X,F9.5,3X,F30. 29) 

END 
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APPENDIX E. PROGRAM FOR CALCULATION OF 


CHARACTERISTIC IMPEDANCE OF SUSPENDED SUBSTRATE 


+ tt 4 


25 


LINE 


THIS IS AN EXAMPLE PROGRAM TO COMPUTE THE DYNAMIC 
CHARACTERISTIC IMPEDANCE OF SUSPENDED SUBSTRATE LINE 
IN SECOND ORDER APPROXIMATION WITH W=1.25 & ER=9.8 








REAL EEFF ,OMG,FREQ,BETA,W,G,ER,T,A,Z0,Z211,Z1,V,R 
PARAMETER ( W=1. 25E-3, G=1.575E-3, ER=9.8, T=0.255E-3 ) 
PARAMETER ( A=3.2E-3 , BETA = 293.61813 , PI=3. 141592654) 
NOTE : PARAMETER A IN FIGURE 1 IS 2A = 3.2E-3 METER 


FREQ = 9. 78036 

V = -1. 7860-2. 2035*(T/G)+0. 4750*(A/G) 
R = 1. 0835+0. 1007*(T/G)-0. 09457*( A/G) 
Z11 = 6*G/W + (14+4*((G/W)**2) )*0. 5 
Z1 = 60*(V+R*LOG( 211) ) 

OMG = FREQ*PI*2. E+9 

EEFF = (BETA*3. E+8/OMG)**2 

ZO = 21/(EEFF**0, 5) 


WRITE( 6,25) Z0 
FORMAT(’ ',5X,F9. 3) 
STOP 

END 





APPENDIX F. SOME MATHEMATICAL CALCULATION 


These equations are used in computer programs 
1. The equations (Ip,q,m) on page 7 

Yj=jwt, LZ, =jou, Kj = 0,8, 

Y, = jwet, Z,=jouu, Kz = oy .u,£08, 
Y;=jwt, L£3=fou, K3 = @u,£, 


Note that Y, = Y; Z,=Z, K? = K} 


. The equations (1n,o) on page 7 


fi Y JOE, ss 1 JO“, 
9 a? +p? — K? V2 a 4p? KS 
y od 7 y => 
ye", JWEE, a 25 JOU, 
2 2 2 2 2 72 
3 a” + B° — K3 es a + p°—K 
Yy3 3 JOE, 723 = 23 JOU, 


. The equation (1g) on page 6 


= NE 
oO DE 


Where : 
NE = Vy2Ct3 + Vy3Cty 
Ja? + p? — Kj 0? + B? — K} 


NE 
JOE,E, JOE, 


Ct 
NE = Le + PK ECR a Ep Ke 


x Ct, 
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(1) 
(2) 
(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 





Ct 
NE=——[ > x yp + Cty x 73] 














JE, 
Defined — YY", then 
JME, 
ee a Sl 
Ee deg 2nfe, 


YE = Complex(0.0, YYE) 
DE = DE\ + DE2 


DE| = ChCt; + CtCr,( 


JWE,E, JWE, 
2 2 2 2 2 2 
Yy3 an + p = K3 ao + B = K3 
ice, 3 
Yyl JWE, JWE,E, 


Since y,3; = 7,1, simplicity (12) and (13) 
\/ a + P= -— K? 
DE\ = CtCts + aes ae Tea 
Je +p? - 
! 


DE\ = ChCty + CHC xb 


af a? + p? ~ K3 
DE£2 = Ct, Ct, —————_— * E, 
a? + p? — kK} 


~ 
2|- 


DE2 = CC xe, 


DE = Ct,Cty + Ct,Cts x = x + + ChCt, x 2 x E, 
In the programs use a symbols: 


COT1, COT2, COT3 for Ct, Ct, Ct,, and 
GM1, GM2, GM3 for 71, 72, ¥3- 


. The equation (1i) on page 6 


Ct, = Coth yh = Coth {,/07 + B? — K? xh} 


IF (a? + B? — K}?) < 0.0, then y, is complex. We defined y, as G,, so 
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(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 





1 
Tanh iG\h 


Ct, = Coth JG,h = 


rea Tal 
Ct, = Complex(0.0, Sac ) 
ELSEIF 

Ct, = Coth Gh = PRIGT 

Ct, = Complex( ei , 0.0) 


The same process for Cr, and Ct, in equation (1j) and (1k). 


. The equation (1h) on page 6 


Where : 
NH = ¥22Ch + ¥23C83 


sar tap RA ja? + p? — K} 
————— i 


NH = ; 
JOU oH, Jou, 


x Ct; 





; 1 Cty 2 no ae 
NH=>,,—L i fo? + B® — K3 + Ctz,/a? + B? — K} J 


We assumed that the value of w, = 1.0, then we have 
1 
NH =——[Ct, x y2 + Ct; x 93] 
jou, a2 3% ¥3 


Write NH in a symbolic form, 


1 
Jou 





N= x [&&& J 


2 
DH = yyV22Cty Cty + 92,9230 Cts + ¥22723CN CL + V2 
Jot PP KP for +f — K3 


DHi = x >—_———_ x CC 
JOU, JOUM, 
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(26) 


(27) 
(28) 


(29) 








— 
ap? Ke Jor + PP — K3 


DH2 = ————_— x CCB 
JOU, JOU, 
<j +f? — KP Jot BP KS 
DH3 = —————. x ————_—_ x C2,Ct; 
JOH, JOU, 


Se 

Ver P- 8 
JOU Ly, 

Simplify (29 to 32): 


DH1 = | ~ x \/0? + p? — Kj x Ja? +p? — KZ x ChCh 


DH4 = ? 

















(ou,) 
DH2 = ; 5 So EPS Ke ea fae eh Ke & CHC 
WL, 
DH3 = a : x fa? + Bp? — K3 x \/ a" + B? — K3 x Ct>Cty 
Ho 
DH4 = a 5 x fa? +p? — K2 x Ja? + p? - K} 
DH = = (DH1 + DH2 + DH3 + DH4) 
; 


Write DH in a symbolic form, 


| 


(jou)? 


DH= x [sss] 





We defined in the program that K? = K? = AKI and K3 = AK2 
aja? + Bp? — KP x Ja? +f? — Kj = AKKI 

a? + B? -- AK, = AKKI 

and 


la? + PP — KZ x ./a? + B? — KE = AKK2 


a? + B? — AK, = AKK2 


From these equations we have the value of ZH as the following: 
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(33) 


(34) 


(35) 


(36) 


(37) 


(38) 


(39) 


(40) 








altos [&&e& | 
JOU, 





[&&&] 


Eanes oe (41) 
(—_}*x Esss] Whe *™ sss] 
jou, 


ZH = 


Defined jou, = YYH 


YH = Complex(0.0, YYH) (43) 
So 
ZH = YH x [&&&] 


[sss] 


75 





— 








LIST OF REFERENCES 


. Tatsuo Itoh, “Anvlysis of Microstrip Resonators”, JEEE Transactions 


on Microwave Theory and Techniques, Vol. MTT-22, November 1974. 


. Tatsuo Itoh, Numerical Techniques for Microwave Passive Structures, 


pp. 334-342, John Wiley & Sons, 1989. 


. Tatsuo Itoh and R.Mittra, “A Technique for Computing Dispersion 


Characteristics of Shielded Microstrip Lines”, JEEE Transaction on 
Microwave Theory and Techniques, October 1974. 


. E. Yamashita and R. Mittra, “Variational Method for the Analysis 


of Microstrip Lines”, IEEE Transaction on Microwave Theory and 
Techniques, Vol. MTT-16, April 1968. 


. E. Yamashita, “Variational Method for the Analysis of Microstrip- 


Like Transmission Lines”, JEEE Transaction on Microwave Theory 
and Techniques, Vol. MTT-16 No.8, August 1968. 


. P. Silvester, “Equivalent Capacitances of Microstrip Open Circuits”, 


IEEE Transaction on Microwave Theory and Techniques, Vol. 
MTT-20, 1972. 


. Harry A. Atwater, Introduction to Microwave Theory, Robert E. 


Krieger Publishing Co., 1981. 


. Rome Air Development Center, “Computation of Lumped Micro- 


strip Capacities by Matrix Methods - Rectangular Section and End 
Effect”, IEEE Transaction on Microwave Theory and Techniques, May 
1971. 


. Yong-hui Shu, Xiao-xia Qi, Yun-yi Wang, “Analysis Equations for 


Shielded Suspended Substrate Microstrip Line and Broadside- 
Coupled Stripline’, JEEE MTT Digest, pp. 693-695, 1907. 


10. Fred E. Gardiol, “Careful MIC Design Prevents Waveguide Modes”, 


Microwaves, pp. 188-191, May 1977. 


11. David Rubin and Alfred R. Hislop, “Design Techniques for Sus- 


pended Substrate and Microstrip”, Microwave Journal, October 1980. 


12. Robert Grover Brown, Robert A. Sharpe, William E. Post, Line, 


Waves, and Antennas, John Wiley & Sons, 1973. 


16 








13. Peter A. Rizzi, Microwave Engineering Passive Circuits, Prentice-Hall, 
Inc. 1988. 


14. Raj Mittra and Tatsuo Itoh, “A New Technique for the Analysis of 
the dispersion Characteristics of Microstrip Lines”, JEEE Transaction 
on Microwave Theory and Techniques, Vol. MTT-19, No. 1, pp. 47-56, 
January 1971. 


15. E. J. Denlinger, “A frequency dependent solution for microstrip 
transmission lines”, JEEE Transaction on Microwave Theory and 
Techniques, Vol. MTT-19, pp. 30-39, January 1971. 


16. T. Itoh, R. Mittra, and R.D. Ward, ” A method for computing edge 
capacitance of finite and semi-infinite microstrip lines”, ZEEE Trans- 
action on Microwave Theory and Techniques, (Short paper), Vol. 
MTT-20, pp. 847-849, December 1972. 


17. Yong-hui Shu, Xiao-xia Qi, Yun-yi Wang, “Closed form equation 
for Shielded Suspended Substrate Line”, Appendix of [ Ref. 9] 


18. Tomoki Uwano, “Accurate Characterization of Microstrip 
Resonator Open End with New Current Expression in Spectral- 
Domain Approach”, JEEE Transaction on Microwave Theory and 
Techniques, Vol. 37, No. 3, pp. 630-633, March 1989. 


19. A.K. Kotronis, “Millimeter Wave Filter Design with Suspended 
Substrate Line”, M.S. Thesis, 1989, Naval Postgraduate School, 
Monterey, CA. 


20. Clifford M. Krowne, “Relationships for Green,s Function Spectral 
Dyadics Involving Anisotropic Imperfect Conductors Imbedded in 
Layered Anisotropic Media”, ZEEE Transaction on Antennas and 
Propagation, Vol. 37, No. 9, pp. 1207.1210, September 1989. 


77 


INITIAL DISTRIBUTION LIST 


Defense Technical Information Center 
Cameron Station 
Alexandria, VA 22304-6145 


Library, Code 52 
Naval Postgraduate School 
Monterey, CA 93943-5002 


Professor Harry A. Atwater 
Naval Postgraduate School, Code 62An 
Monterey, CA 93943 - 5002 


Professor Jeffrey B. Knorr 
Naval Postgraduate School, Code 62Ko 
Monterey, CA 93943 - 5002 


Professor Tatsuo Itoh 

Department of Electrical and Computer Engineering 
The University Of Texas At Austin 

Austin, Texas 78712-1084 


Professor Tomoki Uwano 

Image Technology Research Laboratory 
Matsushita Electric Industrial Co., Ltd. 
1006 Kadoma Osaka. 571 Japan 


Library of the Indonesian Naval Academy 
AKABRI LAUT Morokrembangan 
Surabaya 

Indonesia 


Library of the Naval Institute of Technology 
PUSDIKITAL, KODIKAL, Morokrembangan 
Surabaya 

Indonesia 





No. Copies 




















Ledr. Sudjiwo 

Jin. KRI. BERUANG No. 17/B 
Kompleks TNI-AL, Radio Dalam 
Kebayoran Baru 

Jakarta-Selatan 12140 

Indonesia 


Chief of the Defense Attache 
Embassy of the Republic of Indonesia 
2020 Massachusetts Avenue N.W. 
Washington, D.C. 20036 


. Director of Education of the Indonesian Navy 


Direktorat Pendidikan TNI-AL 
Mabesal, Cilangkap 
Jakarta-Timur 

Indonesia 


Head of Education of the Department of Defense and Security 
KAPUSDIKLAT Departemen HANKAM 

Jin. Pangkalan Jati No. 1 

Jakarta-Selatan 

Indonesia 








